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ON  THE  DETERMINATION  OF  THE  SHAPE  OF  THE  GEOID 
AND  THE  SHAPE  OF  THE  EARTH  FROM  AN  ELLIPSOIDAL 
SURFACE  OF  REFERENCE. 

ABSTRACT: 

By  Introducing  gravity  data  into  geodesy  it  has  been  possible  to  solve  some  of  the 
geodetic  problems  not  only  in  a  purely  geometrical  way  but  also  with  full  considera¬ 
tion  of  the  dynamic  parts  of  the  problem.  The  first  approach  for  the  determination  of  the 
shape  of  the  earth  by  the  aid  of  gravimetric  data  is  based  on  the  famous  formula  of  Stokes. 

According  to  this  formula  it  is  possible  to  determine  the  shape  of  an  equipotential 
surface  if  the  gravity  is  defined  in  all  parts  of  the  surface.  The  method  has  been  extensively 
used  up  to  now  for  a  determination  of  the  so-called  geoid  of  the  earth,  in  spite  of  the 
fact  that  the  formula  is  only  correct  for  a  spherical  surface  and  not  for  a  spheroid  such  as 
our  own  earth.  Another  objection  to  using  this  formula  is  the  fact  that  one  has  to  know 
the  pavlry  values  on  the  so-palled  geoid  and  not  on  the  surface  of  the  earth.  However, 
this  is  in  practice  an  important  limitation  because  measurements  can  normally  not  be  made 
at  die  surface  of  the  geoid  in  other  places  than  on  the  Oceans. 

Evidently  we  have  in  total  two  fundamental  objections  against  using  the  method  of 
Stokes  for  our  geodetic  problems.  During  the  last  1$  yean  there  has  been  extensive 
work  in  this  field  in  order  to  Improve  the  present  technique.  The  most  Important  contribu¬ 
tions  have  been  made  by  Ruslan  scientists  as  for  example  Molodensky  and  Sagreblp.  Moloden- 
sky  showed  in  his  contribution  that  it  was  possible  to  solve  the  integral  equatioo  which  defines 
the  shape  of  the  earth  directly  from  gravity  measurements  at  the  physical  surface  of  the  earth. 
Sagrebin  made  a  similar  solution  of  the  integral  equation  when  there  is  a  change  of  the 
reference  surface  from  a  sphere  to  an  ellipsoid. 

Although  the  two  methods  have  not  hitherto  beat  used  to  any  wider  extent  for  practical 
applications,  it  is  still  evident  that  they  have  given  a  new  view  of  some  of  the  fundamental 
problems  in  geodesy.  For  a  special  study  of  the  ellipsoidal  surface  of  reference,  the  method 
of  Sapebin'  hat  been  chosen  as  the  theoretical  background.  However,  a  study  of  the  deduction 
of  the  sum  of  certain  Legendre  polynomials  used  b>  Sapebin  shows  that  most  of  SapeMns 
final  formulas  are  incorrect.  Therefore  it  was  considered  important  to  resolve  the  problem  by 


t 
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the  aid  of  new  derivations.  This  study  is  based  on  the  same  formulas  as  Sagrebin  used 
for  the  conversion  of  the  Lame  functions  to  Legendre  polynomials. 

For  a  final  computation  of  the  resolvent  equation,  all  functions  have  been  recomputed 
in  order  to  get  a  correct  resolvent.  The  new  expressions  are  somewhat  more  complicated 
than  those  of  Sagrebin  and  an  electronic  computer  was  required  for  the  final  study. 

According  to  the  method  used  by  Sagrebin,  it  is  necessary  to  make  not  less  than 
nine  complete  integrations  over  the  earth  in  order  to  obtain  the  final  value  of  the  geoidal 
height  of  one  point.  This  means  that  the  method  is  too  tedious  for  practical  use.  Another 
way  to  approach  the  problem  is  to  make  use  of  an  iterative  method  such  as  Molodensky 
has  suggested.  There  is  no  objection  to  such  a  method,  except  that  in  most  cases  it  is 
desired  to  obtain  the  final  answer  after  just  one  integration.  It  is  evident  from  a  purely 
abstract  point  of  view  that  repeated  integrations  will  increase  the  accuracy  of  the  computa¬ 
tions. 

Preliminary  studies  made  it  clear,  however,  that  little  increase  of  the  final  accuracy 
is  obtained  by  using  such  a  high  number  of  integrations,  in  the  method  used  here  for  com¬ 
puting  the  resolvent,  all  integration  steps  are  taken  into  consideration,  but  all  this  informa¬ 
tion  is  compiled  in  such  a  way  that  t*o  integrations  are  sufficient. 
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THE  FIGURE  OF  A  SURFACE  DETERMINED 
BY  THE  AID  OF  GRAVITY  DATA. 


There  are  a  few  fundamental  relation.'  which  are  of  utmost  importance  for  any  geodetic 
use  of  gravity  data.  The  Newtonian  potential  of  a  body  is  determined  by  the  function 

W  =  iM-£-  dV 
r 

where 

W  *  the  potential 

p  ■  the  density  of  mass 

V  *  the  volume 

r  =  the  distance  between  the  volume- element  and  the  actual  point. 


The  potentional  W  is  said  to  be  harmonic  if  it  satisfies  the  Laplace  equation: 
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d  x 
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The  potential,  due  to  the  gravitation  effect  of  mass,  is  harmonic  in  all  points 
not  occupied  by  this  mass.  In  a  continuous  massdistribution  the  potential  satisfies  the  so- 
called  equation  of  Poisson: 


A  Vs  »  -  4tip 

For  a  rotating  body  the  potential  is  no  longer  harmonic  as  the  total  potential  is  de 
fined  by  the  function 


.  ( r  i  o  2  2 

VV  *  '  )y)  -y~dV  +  0.5  u  p 

where  p  is  the  distance  of  the  point  to  the  axis  of  rotation. 
The  gradient  is  now  given  by 


grad  W  «  -tSSS-^ 

cos(r ,  x) 

2 

dV  +  u  p 

COS(p.X) 

r 

cos(r,  y) 

cos(p,  y) 

cos(r ,  z) 

cos(p,  z) 

For  a  simple  solution  of  geodetic  problems  it  is  necessary  to  work  with  an  harmonic 
potential- function.  Therefore  the  geodetic  potential  is  normally  replaced  by  an  harmonic 
parameter  which  is  defined  as  the  difference  between  the  geodetic  potential  and  the  po¬ 
tential  of  a  theoretical  earth.  If  these  two  potentials  are  identical  with  respect  to  the  centri¬ 
fugal  part,  then  our  new  parameter,  the  so-called  disturbance  potential,  is  harmonic,  and 
the  well-known  Green  theorem  is  valid 


where 


Vv  , 


up  -TP  =  [(W-U) 


~~  -  — (  ~  -  ■—  ;*1  ds 
dn  r  dn  dn  J 


Wp  =  the  true  potential  at  the  actual  point 

Up  =  the  potential  of  the  theoretical  earth  at  the  actual  point 

n  =  the  normal  of  the  surface 

r  ■  the  distance  between  the  actual  point  and  the  running  point 
S  *  the  surface. 


in  order  to  obtain  an  harmonic  disturbance  potential  it  is  required 
that  the  theoretical  earth  (the  reference-surface)  and  the  actual  earth 
coincide  with  respect  to  their  axes  of  rotation. 

This  is  the  first  important  condition  for  a  solution  of  conventional  type,  in  order  to  fa¬ 
cilitate  the  solution,  another  important  condition  will  be  added  later.  For  a  full  understanding 
of  this  approach  we  note  that  in  case  we  have  a  given  surface,  it  is  possible  to  compute  the 
potentials  directly  by  aid  of  Greens  theorem.  This  problem  is  the  so-called  Neumann  problem. 

If  the  surface  is  not  known,  then  there  is  no  simple  straight-forward  solution.  The  solution 
according  to  Stokes  has,  however,  been  a  useful  approximation. (See  equation  28  "The  Earth 
Form  from  Potentials  and  Gravity".)  It  has  however  to  be  remembered  that  any  solution  of  this 
type  is  based  on  the  approximation  of  performing  the  integration  over  a  sphere  instead  of  over 
the  true  surface  of  the  earth.  Of  much  greater  importance  is  another  limitation  which  is  caused 
by  using  Stokes'  formula.  From  our  formula  26  we  find  that  the  Legendre  polynomial  of  order 
one  is  omitted  from  our  solution.  This  means  that  the  gravity  centers  of  our 
theoretical  earth  and  the  actual  earth  have  to  coincide.  This  is  the  second 
important  condition  for  a  solution  of  conventional  type. 


5' 


We  have  to  note  that  the  solution  from  the  Stokes  formula  Is  only  valid  for  an  equi- 
potential  surface.  It  can  be  proved  that  for  a  number  of  applications  the  solution  gives  a 
useful  approximation  but  when  it  becomes  necessary  to  make  a  detailed  study  of  the  shape 
of  the  earth,  then  Stokes'  formula  is  insufficient. 

The  following  study  is  made  in  order  to  obtain  a  resolvent  which  can  replace  Stokes’ 
formula  for  an  ellipsoidal  reference  surface. 
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STATEMENT  OF  THE  PROBLEM. 


The  following  quantities  are  assumed  to  be  given  at  the  outset: 

a)  Some  reference  surface:  In  the  general  case  this  may  be  an  ellipsoid  with  three 
uaequal  axes,  but  because  of  the  very  small  eccentricity  of  the  equanr  this  will 
later  be  taken  as  an  ellipsoid  of  revolution. 

b)  Gravity  anomalies  Ag  ■  g  -  y  where  g  is  the  acceleration  of  gravity  on  the  geoid 
and  y  is  the  acceleration  of  the  theoretical  gravity  on  the  reference  surface. 

The  quantities  sought  are  the  height- differences  N  between  the  geoid  and  the  reference  sur¬ 
face,  measured  along  die  normal  of  the  reference  surface.  The  height- differences  N  are  related 
to  the  disturbance  potential  T  according  to  the  following  formula: 


N  = —  (1.1) 

r 

The  geoid  can  be  determined  as  soon  as  an  equation  is  obtained  to  solve  the  T.  If  the 
reference  surface  is  a  sphere,  the  required  relation  is  given  by  the  so-called  "fundamental 
equation  of  gravimetry": 


2T  a  T 
R  +  d  n 


which  is  satisfied  by  Stokes' solution: 


-  Ag 


(1.2) 


T  *=  S  Ag  S(t)  dS  (1.3) 

Where  n  is  the  outer  normal  of  the  reference  surface  and  dS  is  the  surface  element  of  a  sphere 
with  radius  R.  The  goal  is  now  to  obtain  a  similar  boundary  condition,  valid  for  an 
ellipsoidal  reference  surface. 
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THE  FUNDAMENTAL  EQUATION. 


In  dealing  with  problems  concerning  ellipsoids,  it  is  advantageous  to  use  the  so-called 
"ellipsoidal  coordinates".  The  orthogonal  system  of  surfaces  is  a  set  of  confocal  qua¬ 
drics,  represented  by  the  equation 


2  2 
X  -  h 


x2-k2 


■=  1 


(2.  >) 


Considering  (2. 1)  as  an  equation  in  X  one  can  distinguish  three  sets  of  roots,  depending 
on  the  sign  of  the  denominators  in  equation  (2. 1 ).  Calling  these  sets  p,/u ,  v  we  have  the, 
relation 

2 


2.2  2,2 

p  ;■>  h  >  ft  >  k  >  w 


(2.  2) 


from  (2.1)  and  (2. 2)  it  follows  that  the  p  defines  a  set  of  ellipsoids.  It  must  be  no¬ 
ticed  that  the  three  axes  of  these  ellipsoids,  in  order  of  increasing  magnitude,  are  located 
respectively  along  the  x,  >'  and  z-axes.  /i  and  v  define  two  sets  of  hyperboloids, 
orthogonal  to  ih,:  ellipsoids  and  to  each  other. 

Equation  (2.1)  must  be  transformed  in  order  to  make  the  quantity  h  equal  to  unity. 


x  ‘ 

1T) 


i±L  HK 


4'' 


X  2  i,  2  X  2 

(T>  -<T>  <T> 


This  can  be  rewritten  as  follows: 


(2.  3) 


x2-. 


+  2  2  + 
x  -  n 


=  i 


(2  4) 


Therefore  relation  (2.2)  becomes 


2  2  2  2 
p  >  '  >  /“  >7  >  v 


(2.  ft) 


Every  point  in  space  is  determined  by  giving  the  ellipsoidal  coordinates  p,/.  ,  v. 
tf  *j>  *2  an**  b  »re  the  half- axes  of  the  reference  ellipsoid  in  question,  in  order  of 
decreasing  magnitude,  this  ellipsoid  can  be  characterized  by  the  relations 


P  =  P, 


(2.6) 


and 
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(2.  6) 


Let  U  be  the  gravity  potential  of  the  reference  ellipsoid  which  at  the  surface 

p  =  p  has  the  value  U 
K  r  o  o 


=  (27) 
p  =  po  O 


The  acceleration  of  gravity  on  this  surface  is  determined  by  the  equation 


.  a  u  _  .  a  U  dp 

y  '  dn  '  'dp  p=f 


(2.8) 


The  gravity  potential  of  the  geoid  is  denoted  by  W.  W  is  considered  to  be  composed  of 
the  potential  of  the  reference  ellipsoid  U,  and  an  additional  part,  the  disturbance  potential  T. 
On  the  surface  of  the  geoid 


W  =  U  +  T  =  W 


(2.  9) 


The  acceleration  of  gravity,  g,  is  obtained  by  the  formula: 


8  .  .  .  5LI  .  .  m  $L  .  *1  (2.1 0, 

6  a  n  3n  dn  dp  dn  dn 

If  Ap  is  the  increment  in  the  ellipsoidal  coordinate  p,  corresponding  to  an  increment  N 
along  the  normal  of  the  reference  surface,  the  following  relation  is  valid 

Ap  *  N  (,2.  11) 

d  n 

In  (2. 10)  we  replace  U  by  the  first  two  terms  of  Taylor's  expansion: 


or 


d  U 

U  s  (u>  ♦  Ap  (-i-) 

P-Pq  9P  P-P  o 


(U) 


P-P  o 


♦ 


N  ( 


a  u  ap 

a p  dn  'p  =  p0 


(2.  12) 


(2. 13) 
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2 

The  next  term  can  be  shown  to  be  of  the  order  of  a  (a  is  the  flattening  of  the  earth). 

2  -5 

Terms  of  this  order  will  be  neglected  in  all  subsequent  formulas,  (a  « 10  ). 

(2.13)  is  reduced  with  the  aid  of  (2.7)  and  (2.8)  to 


U  =  U  -  Ny 
o 


(2. 14) 


Substituting  (2.14)  into  (2.9)  gives 


U  -  Ny  +  T  =  W 
o  o 


(2. 15) 


T  =  W  -  U  +  N  > 
0  o 


(2.16) 


If  we  assume  W  =  U  we  have  found  equation  (1.1).  This  relation  is  sometimes  called 

O  0 

"the  lemma  of  Bruns".  We  now  proceed  to  derive  the  fundamental  equation  for  T.  By  Taylor's 
expansion  of  (2. 10)  we  obtain 


,  3U  .  a2u 

8  - )  ♦  AP(  - ?) 

dP  p«p.  dP  p-po 


X 


dp  d  dp 

( - )  +  AP( - ) 

d n  p * p0  dp  dnp  = 


Po 


a  t 
a  n 


As  in  (2.12)  we  neglect  here  the  terms  with  Ap  so 


.ApJ  (i|^, 

dp  dn  dp  dn  v  dp  dp  dn' 

P  =  P  o  k  P'~P0  P  =  P  o. 


a  t 
a  n 


using  (2.8)  and  (2.11)  we  get 


a»  a t 

g  *  y  =  N  —4  -  — — 
6  '  d  n  3  n 


(2.17) 


and  taking  (1.1)  into  account 


T  ay  dT 

8  '  y  =T^‘  77 


(2.18) 


In  expression  (2.18)  the  derivatives  of  y  and  T  are  both  taken  along  the  normal  to  the 

ellipsoid.  This  is  not  exact  in  as  much  as  T  must  be  differentiated  along  the  normal  to  the 

geoid.  In  neglecting  the  angle  c  between  the  two  normals  (called  plutnbline  deflection)  we 

introduce  only  a  small  error,  namely  the  difference  between  cos  e  and  1.  Since  e  is  always 

3 

smaller  than  one  minute  of  arc  the  error  is  of  the  order  of  a  .  We  can  safely  use  formula 
(2. 18)  or  in  its  final  form: 
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1  dy  „  fl  T 

T  T  +  ~  '  Ag 


(2.19) 


This  important  formula  can  be  derived  from  simple  geometric  reasoning. 

0 


Fig.  (2.1) 


N  | 

_ 

IP 

Ellipsoid  of  reference. 


Omitting  terms  of  the  order  of  N  ,  we  obtain  for  the  theoretical  gravity  at  the 
point  0: 


y  (Q)  =  y  ( P)  +  N 

o  n 
or 

y  (Q)  •  y  (p)  ^  (2.20) 


On  the  other  hand 


g(Q)  = 


a  K 

d  n 


On  / 


a  t 

d  n 


7  (Q  )- 


a  t 
a  n 


which  with  the  aid  of  (2.20)  results  in  equation  (2. 19). 
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DE  T  E !'  M I  NA  T 10  N  OF  THE  UNDULATIONS  OF  THE  GEOID 
WITH  RESPECT  TO  THE  ELLIPSOIDAL  REFERENCE  SURFACE. 


We  now  turn  to  the  solution  of  the  boundary  value  problem,  defined  by  (2.19).  In 
this  solution  we  are  going  to  use  a  rotational  ellipsoid  as  the  reference  surface,  although 
some  of  the  intermediate  formulas  are  referred  to  the  more  general  case  of  a  tri- axial 
ellipsoid.  Hence 

a  =  a  =  a , 

1  2 

According  to  a  formula  of  Bruns  (1)  the  cc efficient  of  T  in  (2. 19)  is 


a_£  ( _L  +  _L>  +  2^2 

y  dn  +  pn  +  y 


(3.  1) 


Herein  is  the  angular  velocity  of  the  earth.  p|n  and  pn  are  the  principal  radii 
of  curvature  at  the  point  of  the  ellipsoid. 


m 


3_ 

2  2  2 
a  (1  +  e  sin  8) 


I  + 


(3.  2) 


a  (I  +  e  t  sin  8  ) 


(3.3) 


2 


b 


2 


where  e}  is  called  the  second  eccentricity  of  the  ellipsoid.  8  is  the  reduced  latitude  of 
the  point  in  question. 

The  boundary  relation  for  T  can  be  written  in  the  form 


(3.4) 


(1)  Bruns,  H.  Die  Figut  der  Erde.  Publik.  Pretisi:.  Geodlt.  Inst.  Berlin  1878. 
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Here 

Ag*=  dp  =  \J  '  +  C)  s  1  n  ^  ®  ‘  ^  8  (3.5) 

4  a  -  b 

Neglecting  ternu  of  the  order  of  ef  and  since  a  - - - -  or 

b  =  ( 1  -  a)  a ,  we  have 


and 


2  a  (2  -  a)  _  .2 

e  ,  =  - - - „  «  =  2  a  + 


!  (1  -  a)2 


0(a  ) 


1  d 


y  3p  p 


(1  +  q  +  a  cos  0  ) 


(3.  6) 


where 


d  = 


2 

w  a 


Hence  the  final  boundary  condition  is 


.  ,  2  „  2  T  d  T  * 

(  1  +  q  4  a  cos  0  )  - ♦  — —  =  -  dg 

p  dp 


(3.7) 


This  formula  is  an  extension  of  the  fundamental  formula  of  gravimetry  for  the  case 
of  a  slightly  flattened  ellipsoid  of  revolution.  As  is  easily  seen,  the  relation  (3.7)  reduces 
to  (1.2)  for  the  case  of  a  sphere  (  q  =  o;  a  =  o;  4fr-  p  • 

Since  the  quantities  of  the  order  of  a2  are  neglected,  we  can  write  for  q  the  value 
it  has  at  the  equator 


q  s  • 


2 

w  a 


Writing  C  for  (1  ♦  q)  and  9(0)  for  cos  0,  which  is  a  function  of  the  position  on  the 
ellipsoid,  the  coefficient  between  parentheses  in  (3.7)  becomes  C  ♦  6(0)  .  T  is  expanded  in 
a  series  according  to  powers  of  a,  for  which,  as  previously  explained,  only  the  first  two  ternu 
are  used 

f  O  1  til 

(3.8) 


T  *  T(0)  ♦  a  T(1) 


Combining  (3.8)  and  (3.7),  and  comparing  the  coefficients  of  equal  powers  of  a, 
we  get  the  two  equations 


2T 


(o) 


dT 


(o) 


dp 


- 


(3.9) 
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and 


2T 


(U 


a  t 


(i) 


dp 


■=  -  e  (s) 


2T 


(o) 


(3. 1 0) 


If  the  quantity  is  found  from  (3.9)  the  right  hand  expression  in  (3.10)  is 
known  and  (3.10)  can  be  solved  in  the  same  manner  as  (3.9).  The  complete  disturbance' 
potential  is  then  obtained  from  (3.8). 

The  height- differences  N  are,  according  to  the  lemma  of  Bruns,  also  split  up  into 
two  terms. 

N  =  N(0)t  aN(,)  (3.11) 

In  all  series  developments  in  the  preceding  formulas,  as  well  as  in  tne  following 

O 

ones,  the  terms  of  the  order  of  magnitude  of  a  and  higher  orders  are  neglected.  The 

o 

error  in  the  height- differences  N  may  thus  be  of  the  order  of  a  N,  or  of  the  order  of 
1  cm.  Errors  of  this  order  can  always  be  neglected.  In  Stokes’  solution  the  flattening  of 
the  earth  is  entirely  neglected. 

The  solution  resulting  from  the  following  method  shows  an  improvement  by  a  factor  of 
a  over  the  classical  Stokes'  approach. 
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EVALUATION  OF  THE  DISTURBANCE  POTENTIAL. 


In  the  following  derivations,  use  is  made  of  the  previously  mentioned  ellipsoidal 
coordinates  p,jt  and  v,  and  of  some  special  harmonic  functions  of  p,  and  v,  denoted 
by  R  (  p),  M  fys)  and  N  ( v),  the  so-called  Lamt'  functions  or  ellipsoidal  harmonics.  (The 
reader  who  is  unfamiliar  with  Lamp’s  functions  may  consult  any  convenient  treatise  about 
this  topic  ) 

The  problem  is  to  find  a  function  T^  which  is  harmonic  outside  of  the  ellipsoid, 
p  =  po>  and  which  on  the  ellipsoid  still  satisfies  the  boundary  condition  (3. 9).  A  bar  is 
added  to  the  quantities  on  the  ellipsoid  to  distinguish  them  from  those  in  space.  The  right 
band  member  of  (3.9)  must  be  regarded  as  a  predetermined  function,  -f  {ft.,  v  ),  of  the 
two  ellipsoidal  coordinates  ft  and  v. 

f  (ft.,  v)  *  Ag*=-^jj-Ag  (4.1) 

d  o 

but  —f  —  can  be  written 
dn 

-J5--  >o V"!  •  "’V  •  w  w 

where,  according  to  Liouville 


2 


2  2,2  2 

*  /i)(p  -  v 


Rj  and  Rg  are  the  first  two  functions  of  Lam£  (see  Poincarl  "Figures  d'equllibre  d'une 
masse  fluide"  1902). 

The  boundary  condition  (3. 9)  becomes 


2T*0)  JT  0 

C -  *  — - — 

p  dp 


f  (ft-,  v) 


where 


f  (M.  v) 


1  R.( p  )  R„(p  ) 
o  1'ro  2  r o 


Ag 


(4.3) 


(4.4) 


The  function  f  (ft ,  u)  can  be  expressed  as  a  sum  of  products  of  the  form  M  (ft;  N  (v) 

as  -j—  is  a  limited  and  continuous  function  on  the  surface.  The  Lame  functions  S,  R,  M  and 
o 

N  will  be  given  two  indices,  n  and  m,  the  first  one  to  denote 
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the  order  of  the  function  and  the  second  for  the  number  within  the  order.  For  the  func¬ 
tions  of  the  first  order  R° ,  R1  ,  R?,  keep  the  notation  used  by  Poincare:  R  ,  R„ 

111  I  4 

and  Rg. 

Although  the  starting  formula  (3.9)  has  been  obtained  for  an  ellipsoid  of  revolution, 
in  this  section  the  general  notation  defining  Lame's  functions  for  a  tri-axial  ellipsoid 
will  be  used.  It  can  be  noticed  that  for  an  ellipsoid  of  revolution  R  =  R  . 

*•  ^  i  \ 

(o) 

Since  T  satisfies  the  equation  of  Laplace  outside  the  ellipsoid,  we  can  express  T 
as  follows: 


t 

n=  2 


2n 

Z 

m  =  o 


A m  Rm(p  ) 
n  n  r  o 


Sm(p)  Mm  (ft)  Nm(i/) 
n  n  ^  n 


(4.  5) 


The  disturbance  potential  T  is  not  uniquely  defined  without  knowledge  about  the 

position  of  the  reference  ellipsoid  and  its  relation  to  the  earth.  Therefore  we  add  the 

assumption  U  =  W  ,  which  implies  that  the  term  of  order  n=  o  is  equal  to  zero, 

o  o 

and  the  assumption  that  the  centers  of  gravity  of  the  earth  and  the  reference  ellipsoid 

coincide,  which  implies  that  also  the  term  of  order  n  =  1  is  equal  to  zero.  As  can  be 

seen  in  (4.5)  the  summation  over  n  is  taken  from  two  to  infinity.  The  function  fi/i.  v) 
is  thought  to  be  expressed  in  the  form 


f  (ft.  v) 


oo  2n 
=  L  Z 
n*2  m*o 


Bm  Mm 
n  n 


(/*■)  N 


m 


(") 


(4.6) 


For  the  sake  of  brevity  the  following  notation  is  used. 


m  m  . ,  m  ,  _  m 

A  M  (^.)  N  (v)  =  P 

n  n  '  n  n 

_  m  „.m  ,  ,  „,m  ,  x  sliti 

B  „  M  (ft.)  N  ( v )  =  g* 

n  n  n  n 

Expressions  (4.5)  and  (4.6)  become  now 


T<o) 


oo 

Z 

n»2 


2n 

Z 

m  =  o 


_  m  ,  ,  „  m  ,  _m 

Rn  Sn  (p)  Pn 

non  n 


(4.7) 


and 
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m  2  n 

-k.  1TI 

*8  =  2  2  S*„ 

n=  2  m  =  o 


(4.  8) 


At  the  surface  of  the  ellipsoid  p  =  pQ  we  have 


H  !<«)."  e"  —  sm  (p  )  r”  (p  )  p" 
p  „  p  n  ro  no  n 

ro  n=  2  m  =  o  ro 


(4.  9) 


=(o)  w  2n 

—  =  £  £  R  (  p  ) 

dp  _  n  ro 

r  n=2  m  =  o 


(^} 


dp  /P=P0  n 


(4. 10) 


From  (4.3).  (4.8),  (4.9)  and  (4.10)  we  obtain 


oo  2n 

E  E  Rm  (p  )  Sm  (p  ) 
„  n  o  n  ro 

n=2  m=o 


L». 

po  sm(P  )  \  aP  frmi) ol  n 

n  ro 


OO  2n 

_  _  -*;n 

°  '  2  2  »„ 

n=2  m=o 


(4.  11) 


Thus, 


a7* 


12C  ♦  1 

/•*n(P>\  1 

V°  s>°> 

\  Op  J  f 
P  =  P0J 

(4.12) 


Substituting  (4. 12)  into  (4.7)  we  get  for  the  disturbance  potential  the  following  ex: 
press  ion 


=  E  Z 


2n 
E 

n*2  m=o 


. m .  ,  *m 
sn  <P)  8 n 


f 2C  1  /a  s  p  (P)\  1 

|/o  +  Sm(p  )  \  /  j 

l.  n  o  dp  p  =  p  / 


-  m  .  , 

S  (p  ) 
n  ro 


(4. 13) 


Replacing  R  (pQ)  and  S  (pQ)  by  R  and  S  we  get  for  P*f>0 


t<o)  »  r 


2n 

E 


5  n 


n  =  2  m*o  l  d  S  _  2C 
- : _ _ _ B - 


-n- 


T<0)  -  p 


oo  2  n 

£  £ 
n-2  m  =  o 


m 


dlnS 

_ n 

d  In  p 


-  2C 


In  the  special  case  of  a  nonrtotating  sphere  this  yields  Stokes' solution. 


(4.U) 


18- 


THE  HEIGHT- DIFFERENCE  N  EXPRESSED  AS  AN  INTEGRAL. 


We  try  to  get  now  an  expression  for  the  height  difference  N  analogous  to  the 
expression  of  Stokes.  Taking  into  account  the  orthogonality  of  the  LaniJ  functions, 
we  can  write: 


g*m  •  — 1  Ag*  Tm1" 
°  4*  Im  '  °  L  " 


where 


m  (,s)  N®  (u)  M™  Os*)  N™(u*) 


d  w  (5.1) 


Tm  ’Cl  J 

I  =  -  )  1  ( M  N  )  d  w 

n  4^  o  n  n 


(5.  2) 


dui  represents  here  the  surface  element  of  the  ellipsoid.  The  series  (4.14)  is  written 
in  a  shorter  form  as 

c«a  2n 


T^°^=p  E  E  T 
n*2  m  =  o 


m 


(5.3) 


,  m 


I? 


j  •••  ^ 

n  Em 


(5.4) 


Thus, 


.  m 


d  in  S  , 


d  In  p 


-  2C 


(5.5) 


_m 
T  * 


*  4*  lm  Ero 

n  n 


$  lo  Ag*M®  <M>  n“(i»)  m”(#s')  n”(i/')  d u  (5.6) 


Using  (6.3)  we  have 


*(0’ .  £  J  , 


OO 

2n 

E 

E 

n*2 

m  -  o 

kJDi  Min  .  ..in  .  t 

M„(^)  N0(,)Mn(M')  n  n( i>) 


IBIB 
n  n 


du  (5.7) 


By  (4.4)  we  obtain 


>g*  = 


1  R,  R„ 

o  1  2 


Noting  that  R  =  p,  we  have 
3 


„(°) _ 


4*  R1R2 


\  igG  (p,#,/,  v)  dw 


(! 


where 


00  2n  mV)n^)mb(;)n>') 

,  ,  n  n  n  n 

Gift,  v  ,fi,  v)  =  i  £  - 


( < 


n  =  2  m  =  o 


jin  Em 

n  n 


For  die  case  of  an  ellipsoidal  reference  surface  we  have  a  formula  analogous 
to  that  of  Stokes: 


>1 


(o) 


4*  rir2  r 


S  Ag  G  {M  ,v,u>\  u)  dw 


or 


\  1  Ag*G  (h.v  m.v)  dw 
4*y  0 


We  have  the  relation  1  dw  >  do  ,  where  do  is  the  surface  element  of  the 
o 

unit  sphere,  and  dw  is  the  element  of  surface  on  the  ellipsoid. 


•  8) 


•  9) 


5.  10) 


5.11) 


THE  CASE  OF  AN  ELLIPSOID  OF  REVOLUTION. 


In  the  courte  of  deiiving  the  expression  for  the  height  difference  N,  the  deviations 
of  the  geoid  from  ellipsoidal  shape,  we  have  not  required  the  ellipsoid  to  possess  any 
rotational  symmetry.  As  pointed  out  before,  the  eccentricity  of  the  equator  of  the  earth 
is  not  very  pronounced.  As  a  reference  figure,  which  docs  not  necessarily  claim  to  be 
the  final  solution,  an  ellipsoid  of  revolution  serves  excellently. 

In  this  case,  the  products  M  N  degenerate  into  ordinary  spherical  harmonics,  and 
the  functions  R  and  S  are  transformed  into  Legendre  polynomials  of  the  first  and  second 
kind  respectively  and  of  Imaginary  argument. 

Thus  we  have  for  the  case  of  rotational  symmetry: 

R .  =  R  ;  M=  sin  6 

2  3 

M m  (fi)  =  X m  (cos  0) 
n  ^  n 


k,m  , 
NnOO 


{sin 

COS 


m  X 
os  m  X 


0  =  90  -  fl.  where  S  is  the  reduced  latitude. 


^1  -  p2  =  i  s  :  p  =  \j  1  +  s2;  s  *  ^p2  -  1  ; 

■:  o«>  *  0is>  sT  <is>  =  o“(i  o 

n  n  n  n 


Thus  we  get 


J.  x  _d.  _ds  Ifl  ♦  s2  _d_ 
dp  ds  dp  s  ds 


_ L  (  Dm  ,  -JCOS  m 

'n  -  „  <C0‘  e>{,in  m  x}_ 


do  = 


( n+m  )1  1 


(2n+ l )(n- m)!  c 


(c  =  1  for  m  =  o;  c  *  2  for  m#o) 
m  m 


(6 
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E 


m 

n 


1  ♦  S2  _d_ 

s  ds 


In 


0 


(is) 


2C 


(6.2) 


C  =  1  +  q 


That  is 


G  v)  =  Z  (  e  .  .ft ,  X) 


(6.3) 


where  0^,  X^  now  denote  the  fixed  point  and  0,  X  the  running  point,  and 


Z  (  6  X  .6,  X) 
0  0 


oo  r* 

£  £ 
n  -  2  m  =  o 


,  m 


(2n+  1  )(  n  -  m ) !  c  P  (cosO)costn(X-X  ) 


m  n 


( n+  m )! 


(6.  4) 

Note:  we  sum  over  from  0  to  n  (not  2n)  because  the  rotational  symmetry  introduces 
a  degeneration.  For  a  given  n.  there  now  exist  only(n+1)  linearly  independent  functions. 

Quite  naturally,  the  product  mJJ1  N”1  must  differ  from  the  elementary  spherical  harmonic 

Pn'  (cos  0  )/C0>  111  >  b\  a  constan:  factor,  but  this  factor  is  cancelled.  3s  its  square 

n  jsin  in  XJ 

occurs  both  in  the  numerator  and  the  denominator  of  the  formula  (  5.  9)  for 
G  (>i.  v  ,fi',  v)  since  I  is  in  the  denominator. 

Thus  (5.  11)  takes  the  form 

N<0)  ’-^7  S  AgZ(VV  e-X)  d«  (6.5) 

^o 

where  Z  it  given  by  (6.4). 

Now  we  transform  the  domain  of  integration  from  that  of  the  ellipsoidal  surface  to 
that  of  the  unit  sphere.  Remembering  that  d  or  =  do  and  taking  (5.11)  into  accouut, 
we  obtain  the  formula: 


22 


4* 


-  $  Ag*Z  (»o.  xo.  e.  X.)  d< 


(6.6) 


#  1/ 2 

When  finally  we  Introduce  the  facts  that  Ag  *  U  l+e1 
we  obtain  for  the  ellipsoidal  surface  of  reference: 


2  2 

sin  6  Ag  and  a 


n(0)  =  777  i  RT^AgZ<VX,  e.  X)  do 


(6.1) 


where  Z  is  given  by  (6.4),  Er  a  determined  by  (6.2).  It  can  be  shown  that  (6.6) 
contains  Stokes'  formula  as  a  special  case,  and  thus  (6. 6)  can  be  regarded  as  an  extension 
of  tha  formula. 

The  function  Z  given  by  (6.4)  must  now  be  studied  more  closely.  Doing  tha  in 
accordance  with  Sagrebln,  we  find  the  following  expression  after  a  series  expansion: 


Z(e0.x0.  6.  X)  =  S(*)  ♦  ef  0  (  *;»  .  x  .  e.  x)  +  2q  *) 


(6.8) 


where  0  and  fr  stand  for  the  following  expressions  : 


OO  _  2 


0(*:  e  ,x  ,e,  x)  =  -  & 


n  ( 2n+ 1 ) 


o  o 


n**2  (n- 1)  (2n+3) 


Pn  (cos  *)  ♦ 


+  2  Z 


2n+1  „  2n+l 

- - -  P  (cos  *>  *  z  - - —  P  (cost) 

n=2  (n-l)*(2n+3)  n»2  (n- 1 )2(2n+3)  6X2  “ 


(6.9) 


t  (*) 


2  p  COf(tj 

n=2 ( n- 1 )®  n 


(6.10) 
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INVESTIGATIONS  OF  THE  FUNCTIONS  0  AND  . 


We  begin  studying  the  function 


eo-  Xo1  0'  =  '  +  2$2  '  #3(*  :  0o'X'o'°-  *■>  (7.1) 

The  first  two  terms  of  (7.1)  are  functions  of  the  angular  distance  *  between  the 
two  points  (  8q.  Xq)  and  p ,  X)  where  6  is  the  complement  of  the  reduced  latitude. 

These  functions  are  expressed  by  the  equations  (cf  6.9) 


<*) 


n2(2n+l) 

“  o 

u=2  ( n-  1 )  *(  2n*3) 


P  (cos 
n 


*) 


(7.2) 


*  00  2n+ 1 

$2<*)  =  £  -  P  (cos  *) 

n  =  2  ( n-  1  )2(  2n+3)  “ 


(7.3) 


Hie  third  term  in  (7.1)  depends  not  only  upon  the  relative  position  of  the  two  points 

but  also  upon  their  absolute  position  on  the  ellipsoid.  Thus  the  third  term  is  a  function  of 

all  four  coordinates  6  ,X  ,9,  X. 

o  o 

If  we  change  the  order  of  summation  and  differentiation  in  the  last  term  of  (6. 9; 
we  have  the  formula 


#3<*ieo-*o,0,X>  “  - 2  *2  <*> 

d  X 

o 


(7.4) 


,  •*♦*(*)  a*  2  «$„(♦)  a2* 

(*!  6o'X  o’  e’  K)  “  -  ( - )  ♦  - - 

3  0  0  a  *2  ax  a  *  ax* 

0  T  o 


(7.5) 


d  t 

It  is  clear  that  the  derivatives  —  and  — —  must  be  functions  of  both  sets  of 


coordinates  0  ,  X  ,  6,  X. 
o  o 


ax 
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The  derivatives  can  be  obtained  from  the  formula  (the  spherical  cosine  theorem): 


cos  ♦  =  cos  0  cos  0  +  sin  6  sin  6  cos  (X-X  ) 
ooo 


(7.G) 


Htus  for  the  evaluation  of  the  function  0  it  is  sufficient  to  obtain  the  functions 
0.  (  *)  and  0O  (+)  whereupon  0  can  be  found  by  ordinary  differentiation.  To  obtain 

I  c  3 

0  and  0n,  we  decompose  the  fractions  preceding  P  (  cos  + )  in  (7.2)  and 
1  £  n 

(7.3)  into  partial  fractions.  We  obtain 


n  ( 2n+ 1 )  _34  J _  +  _3_  1  _  J_8_  _ 1_ 


(n-  1)2(2n+3) 


25  n-1  5  ( n-  1  ) 


25  2n  +  3 


and 


2n«1 


(n- 1  )^(2n+3) 


4 

25  n-  1 


1 _ +  _3_ 


5  (n-1)4 


25  2n+3 


Ulus 


®  26  1 

0  ( +)  +  2  0  (  *)  =  -  E  P  (cos  *)  -  -  Z  -  P  (cos  *) 

’  n=2  "  25  n=2  n-1  " 


3  ® 

+ -  Z 


00 


5  n=  2  ( n- 1 ) 


2  n 


P  (cos  *)  + 


25  n=2  20+3 


P  (cos  *) 
n 


The  function  0(*)  (cf  6.10)  takes  the  form 


oo 

0(t)  =  2  Z 


CD 


P  (cos  ♦)  +  3  E  — -  P  (cos  *) 

n  =  2  n-1  "  n=2  (n-1)2  " 


Writing  the  above  expressions  in  a  shorter  form,  we  have 


26 


0,  <♦)  ♦  2  02(*)  -  -  *0<*)  ‘  I? 


(7.6) 


V*)  *,<♦>  'A  *3  {t) 


(7.  7) 
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♦1 

(*> 

+  3  «>2  (*) 

(7. 

8) 

nee 

we  have  introduced  the  expressions: 

(*) 

=  £ 

P 

a 

(7. 

9) 

(*) 

=  £ 

-~r  P 

(7. 

10) 

n-  1  n 

(*) 

=  £ 

1 

( n- 1 )®  Pn 

(7. 

11) 

(*) 

=  £ 

-L-  P 

2n+3  n 

(7. 

12) 

AU  summations  are  taken  from  n  =  2  to  n= 00 ,  and  Pq  stands  for  Pn  (cos  *). 

The  sums  (7.9  -  12)  can  be  obtained  by  utilizing  the  generating  function  for  the  Legendre 
polynomials: 

1 


00 

E  Pn  (cos  *  )  x‘ 
n*o 


(1  -  2x  cos  ♦  +  x‘ 


or 


00 

£ 

n=2 


1 


£  '■* 


-  1  -  X  cos  * 


1  -  2x  cos  *  ♦  X 
Putting  x  =  1  in  (7. 13),  we  have 


(7. 13) 


1 


f 


1  t 

1  -  cos  *  * - cosec  — -1  -  cos  t 


2  cos  t 


The  sums  ^  (*)  and  (*)  we  of  the  general  form:  £  Pn 

We  introduce  the  expression 


_  1  _  n-k 

C  ~k  Pnx  “  Sk  *> 


By  differentiation  we  find  the  relation 


(7. 14) 


(7.15) 


0  „  , .  ,  _  _  n- k- 1 

r —  S.  (♦,  x)  «  E  P  x 
dx  k  n 
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That  is 


Sk  (*.  x)  =  ^  (  Z  Pnxn‘ k‘ 1 ) 
But  we  have,  by  (7. 13) 


dx 


(1 • 16) 


E  P  x 
n 


n-k-1  1 


k  +  l 


/l  -  2x 


-  1  -  X  cos  ♦ 


COS  *  +  X 


Thus 


Su  (♦•x>  =  S  -V. 


k+1  .  .  . 

o  x  Lr  (x) 


-  1  -  x  ui 


dx 


(7.17) 


(7.18) 


where  r(x)  =  ]/ 1  -  2x  u  +  x2  and  w  =  cos  *  to  achieve  shortness  of  notation 

To  obtain  <p ;  we  put  k  =  1: 


(♦•  *)  =  \  — o( -  -  1  -  X")  dx 

«  x  K*) 


(7.19) 


and  then 


(*)  =  s1  (♦.  1) 

Sj  (*,x)  is  evaluated  as  follows: 


(7.20) 


S,(*.x)  (-= - 

_  x  r(x) 


(Singularities  which  arise  when 


1  w  r  1  1  x 

-  — —  -  — )  dx  =  )  — —  dx  +( - w  In  x) 

x2  x  0  x  r(x)  •  x  o 

x-*o  must  cancel  in  the  end  result.) 


\  — 4^ —  =  (substituting  y  ■  •—)  =  J—l y 
x  r(x)  r(y) 


y  -  u 

r(y) 


but 


* 
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Thus 


y  dy 
r(y) 


[  3  r  ( v )  u  i 

^  ( _ L  +  _  )  dy  =  -r(y)  -  w  In  y  -  u  +  r(y) 

3y  r(y) 


r(x) 


-  u  in  l-xw  +  r(x)  +  w  In  x 


It  follows  then  that: 


S1  (i'*  x)  ~  xr("?-  '  w  In 


1  -  ux  +  r(x) 


K 


As  x  approaches  zero,  - — —■*  ^  approaches  the  constant  w. 


Thus 


S,(+.  X)  .  »  in  ^ 

1  x  2 


1  -  cjx  +  r(x) 


(7.21) 


Putting  x  =  1  we  obtain 


?,(*)  -  .  S  1  (  * ,  1 )  =  1  -  ^2-2  uj  -  ui-  u  In  —  (  l  -  w+1^2  -  2u) 

_  ,  1  -  2  * 

But  ui  =  cos  ♦  ,  thtis  — -  =  sin  — — 

which  yields 

(*)  =  1-2  sin  -  cos  +  -  cos  *  ln(sin2-^-+  sin-y)  (7.22) 

To  obtain  ^(t),  we  must  study  the  more  general  sum 

H  (♦.  x)  =  E  - P  x"'1  (7.23) 

(«-l) 

Differentiating,  we  obtain 


a 

a  x 


H  (t.x) 


l 


I 

n-  1 


P  * 


n 


n-2 


(♦.x) 


(7.24) 


by  the  use  of  (7.15). 
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Thus 


A 

H  (*■  x)  =  $  ~  s,  (*.  X)  dx 


1) 


That  is 


H  (*.  x)  =  $  ( 


Putting  x  =  1, 


JiL  in  — 
x2  x  2 


l 


1  -  OJ  X  +  t(  x) 


-)  dx  (7.25) 


h  (...)  ■ ,,  i.)  •'  (iai>  •  js  - -s-in  i- 

2  0  X2  x  x  2 


We  consider  the  expression 
1 


Ju  (*)  "  \  —  S,  (*.  x)  dx  =  jj^*  u 


In  £  -  w  J 


t  * 


(2) 

le 


where  e  is  a  small  positive  quantity  which  approaches  zero  ,  and 

c  ■  i  ^ 


dx 


e  x- 
1 


dx 

»«  .2 


1  /t 


=  -  1  ♦  1/2-  2w  ♦  w+0(c)-$  ■  -  *  -  ♦  w  $  — ^ 
1  '  *<*)  ^  *(! 


(y) 


1-wx  +  r(x)j)dx 

(7.  2ti) 


(7. 27) 


(7.28) 


*  UJX  ♦ 

I(x) 

I- 

(7.29) 

of  lV> 
l  e 

and 

1 

-5- 

x  -  w 

-  dx  = 

*  J* 

< 

x  r(x) 

I)  Sagrebin  did  not  include  the,  "x"  in  the  denominator  of  the  integrand. 
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1 

where  we  have  used  partial  integration,  the  substitution  y  =  —  ,  and  the  fact  that 


1  -t(Q 
e 


ui  +  0  ( €) 


Remembering  that  )  — — ■  =  In  x  -  w  +  r(x)  +  a  constant  and 
r  (  x )  l  i 

disregarding  terms  O  (e)  (of  the  order  of  e ),  we  have 


=  -  1  +w  +  \j~2~ 


_  ,  1  -  ui  +  Ir2  -  2u7 

2w  -  In  - : - - -  + 

1  -  u ) 


+  ui  In 


-  u  +  ^2  -  &  ui 


-  ui  In  e  =  -1+cos  *  ♦  2  sin  —  - 


l|f  * 

-  (  1  +  cos  +)  In  (1  +  sin  — )  +  (1  -  cos  ♦)  In  sin  —  -  ui  In  < 


Thus  letting  e  approach  zero: 

t  * 

^2(*)  *  -  1  +  cos  +  +  2  sin  —  (1-cos  +)  In  sin  —  - 

*  (2) 

-  (1  +  cos  *)  In  (1  +  sin  — )  -  cos  +  Jv  ' 

&  I 

It  remains  to  evaluate  .  It  is  not  possible  to  express  this  function  in  terms  of 
elementary  functions,  but  it  can  be  transformed  by  the  substitution  x  =  -y-  and  by 
partial  integration: 

1 

ln  |*5*  (>  *  ♦  r(x))j  dx 1 

o 

{  x  =  J_;  ,  = 

y  *  y 

GO 

=  S  J  (In  |  Y  (y  -  W  ♦  r(y)  J  -  In  y  )  dy» 

to 

=  ln2y  ♦  In  y  In  -j-<y  '  «  ♦  *  (y)  )  J 

co 

*  \  •  -7~. )  In  y  dy 

;  y  f(y)  7  7 


.  30) 
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* 


Thus 


L  (+) 


(7 


where  L  (♦)  is  given  by  the  equation: 


L  (*) 


-  $ 


( —7—  -  — )  In  x  dx 
'  r(  x )  x 


(7 


The  term  -°x  -  is  included  to  make  the  integral  converge  at  the  upper  limit.  Thus  the 
expression  for  <p  ^  ( '<')  now  reads: 


♦  + 

<P  =  -  1  +  cos  +  +  2  sin  —  +  ( 1  -  cos  *)  In  sin  —  - 


-  (1  +  cos  ♦)  In  (1  +  sin  — )  +  cos  *  L  ( +) 


(7 


We  now  make  a  closer  study  of  the  function  L  (+).  In  the  limits  for  m  =  +  1  ,  that 
is.  *  =  o°  or  ISO0,  we  may  compute  the  exact  values  of  L  (*)  as  follows: 


♦  =  0  ;  u  =  1 

,(2), 


(o)  =  -  J(,  '(0)  =  -  S 


In  ( 1  - x ) 


dx 


By  series  expansion  we  have 
I 


-o  x»°  *• 

L(o)  =  )(£—)«»*=  t 
9  n=l  n»1 


o 


n- 1 


dx  = 


=  Z 


I  * 


n=  I  n 


.31) 


.  32) 


.  33) 


6 
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2. 


*  =  180  , 

,o 


1 


.1 


L  (.SO’)  -  -  (  d, 

0  * 

L(o)  +  L  (180°)  =  -  $  1  :  JLi 

O  V 


dx 


Of 


Making  the  substitution  x  s  t  ;  we  have  2  In  x  =  In  t  and 
dx  _  1  dt 

x  2  t 


Thus 


L  (o)  +  L(180°)  -i-  L(  o ) 


O  1  n 

L(  1 80  )  =  -  —  L(o)  =  -  — 


The  derivatives  of  L(  *)  can  be  obtained  explicitly  by  differentiating  the  terms 
under  the  integral  sign  with  respect  to  *  (see  7. 32): 


d  L  (  ♦)  d  L(  »)  d  ui 
d  *  "  du  d  ♦ 


-  sin  *  L 


(7.34) 


too 

•  dL  (  x  In  x  . 

L  =  ——  *  )  — -  dx 

dw  i  3,  , 

1  r  (x) 


(see  7.32) 


(7. 33) 


But 


thus, 


d  1 


*  ♦  W 


flx  r(x)  r(x)3  r(x)3 


9  1 


r(x)3  r(x)3  dx  r(x) 


r(x)3  r(x)3 


2  2 
r(x)  2 uix  -  x 


Kx)‘ 


J _  +  otx  x(<u-x) 

.  3  +  . .  ,  3 


r(x)  t(x)  r(x) 


i  us  a  i 

-  +  - o  ♦  x - 


r(x)  r(x) 


dx  r(x) 
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This  gives 

OS 

i  In  x 
L  =  hi  )  - 


dx  ♦  u/2  S 


x  In  x 


1  r(x) 


dx  +  u  y  x  In  x  — — —  dx 


>  Kx) 


a  i 

dx  i(  x ) 


i  .  a  i 

-  )  In  x  -  -  dx 

1  d  x  r(  x ) 


09  i—  -1 90  t 

I  hi  x  ,  2  ,  >  x  1  n  x  i 

L  =  <j  )  -  dx  +  w  L  +  <j  -  -  ui  ) 

1  r(x)  L  r(  x )  1 

09  t—  —  O© 

-  w  i  ii.  .  fULJL  I  *  1  il_ 

1  :(x)  L  r( x )  j,  1  r(x)x 

2  /  r 

=  uj  L  ♦  to  ¬ 


la  x. 


t(x) 


dx 


x  In  x 


i(x) 


-  w+r( x ) |  J  ^  -  j~l n  |—  -  w  +  —  r(x)|j 


Thus 


2  «  1 
( 1  -  w  )  L  =  win  — 


-  w  +  |(2-  2w  1  -  w  +  1/ 2  -  2w* 

- : - +  in  - ‘ - 

2  1  -  w 


2  2  2 

w  In  (s  +  s  )  +  In  (s  +  s  )  -  In  s  =  ( 1  +  w)  In  ( 1  ♦  s) 


(1  -  u)  In  s 


where 


s  =  sin  — 


We  have  now 


'  1  1 
L  =  -  In  (1  +  s)  -  -  In  s 

1  -  w  1  ♦  ui 

and 

d  L  _  /  s  c 

-  =  -  2  s  c  L  =  - In  s  -  —  In  (1  +  s) 

d  *  c  s 

where 

t 

c  -  cos  — 
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Thus  the  formula  for  d  —  ^  is: 

d* 


d  L(  +)  ^  ^  ^  ^ 

-  =  t  g  —  in  sin  —  -  cot  —  In  (1  +  sin  — )  (7.36) 

d  ¥  2  2  2  2 

Thif  formula  is  differentiated  once  again: 


d2L<*)  ,  1  +1  2  *  *  1  2  + 

-  =  1  -  —  cosec  —  +  —  cosec  —  In  (  1  +  s  i  n  — )  +  —  sec  -  • 

d*2  2  2  2  2  2  2 


•  In  sin 


r  .  :<7) 


We  now  turn  to  the  evaluation  of  the  function  (*)  given  by  (7.12). 

O 


V*>  “  1 


2n+3 


1  1 

P  =  —  E - P 

n  n  +  -a-  n 


Putting  k  *  -  —  and  x  =  1  in  the  general  formula  (7.18)  we  have 


1 


1 


*3<*>  =  —  s_3  <*.  I)  “  -M  — ( 

2  -  2  o  K*) 

*•  x 


uix  -  1 )  dx 


That  is, 


1 


*,(  *)  =  — ■  S  p— 

3  -  r 


fx  dx 


2  0  yi  -  2xw  +  x‘  5 


(7.38) 


Thus  it  is  necessary  to  compute  the  integral 

_L  (  _ f*  d* _ 

’  -  1  ]j,  .  * 


2  Q  |/1  -  2xu;  ♦  x 


Substituting  fT *  tg  we  obtain 


(7.  ••i>) 


!,(*)  —  J 

“  2  0 


2  ‘S2fd* 


cos 


J  ^  4  ^ 

1  -  2  w  tg  —  +  tg  — 
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from  which  we  have 


U*>  =  —  ^ 
2  „ 


2  if 

tg  j  d  <p 


,  b2  •  2 

1  -  k  sin  if 


where 


k2  =  ~  ( 1  +  ui)  -  y(1  +  cos  ♦)  or  k 


■  cos  • 


However  this  last  integral  can  be  expressed  as  follows: 


(7.40) 


JL  2  if 
2_.  2  13  2  dv 
o  ^ 


I  i 


2  2 

1  -  k  sin  v 


,  ^  2 

1  -  k  sin  if 


0  ^ 


2  „  t/  ,  .  2  .  2 

0  p-k  Sin  «> 

which  can  he  expressed  in  snorter  notation  as 


\  ^1  -  k2sin2*>  d*> 


J2(t)  =  sin  -y  ♦  -i-  K  (*)  -  E(  ♦) 

* 

where  K  (4}  and  E  (+)  are  complete  elliptic  integrals  of  modulus  k  =  cos  — 
Thus,  finally 


(7.41) 


♦)  *  ‘inT  +  TK(#)  -  E(*)  -  jeos  *  --j  (7.42) 

As  we  are  going  to  require  the  first  and  second  derivatives  of  ?  (x)  in  the  following, 

3 

we  sute  here  the  corresponding  derivatives  of  K  (*)  and  E(t): 


dE 

dk 


K) 


(7.43) 


—  *  —  (  E 
dk  k  1  -  k2 


K  ) 


(7.44) 


Formulas(  7.43)  and  (7.44)  can  be  found  in  many  standard  tables  (  see  Dwight: 
"Tables  of  Integrals  and  other  Mathematical  Eata".  4th  Ed,"  1961,  Macmillan  Company, 
New  York.) 
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Now 


Thus. 


dE  _dk  d_E _ j_  .  *  dE 

d*  d+dli  2Sn2  dk 


dE 

d  't' 


E) 


and  similarly 

=—  tg — —  K  -  cosec  t-E 

d  t  2  b  2 

Differentiating  once  more,  we  obtain  in  the  same  way: 


and 


1 

4 


2  +  „ 
sec  - K 


1  2  * 

-  —  ( 1  +  sec  — )  E 

4  2 


d2K  1  2  #  cos  + 

- -  =  —  tg  —  K  +  - — 

d+  4  2  sin 


(7.45) 


(7.46) 


(7. 47) 


(7.48) 
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EXPLICIT  EXPRESSION  FOR  THE  FUNCTION  Z  (8  ,X  ,  0,X). 

0  0 


The  dominant  function  in  the  exptettion  for  Z  U  the  well-known  function  of  Stokes: 

^  t  t 

S  ( *)  =  00560  —  +  1  -  6  sin-  -  5  cos  *  -  3  cos  +  In  (sin-^- + 

2  2  2 

2  % 

+  sin  -)  (8.  1) 

Considering  the  equations  (6.8)  and  (7.1)  we  can  express  Z  as  follows: 


Z  (8o.  \o.  0 ,  X)  =  S  ( ♦)  ♦  e*  f(*)  -  e®  '0T  (*;  e#.Xo.0.  X)  + 

+  2  q  tfr  (  ♦)  (8.2) 

where  we  have 

«*)  ■-#,<*>♦*♦  2  <*>  (8-3) 

The  explicit  expressions  for  f(  t)  and  ^  (*)  can  be  found  immediately  by  the 
aid  of  formulas:  7.6.  7  8,  7.14,  7.22,  7.33.  7.42. 


♦  (  55  co«  *  — |- )  In  <  1  ♦  sin  -y)  ♦ 

♦  cos  *  ♦  -y)  In  sin  y  ♦  -|-  cos  *  L(*)  ♦ 

♦  8t  k  (t)  •  TS  E  (t)  (8-4) 

•  -  1  ♦  2  sin  —  t  cos  *  -  (3  ♦  5  cos  t)ln(l  ♦  sin -y) 


+  (3  -  5  cos  t)  In  sin -j-+  3  cos  *L  (♦)  (8.8) 

To  obtain  the  function  A  ,  we  need  A  (♦),  which  can  also  be  found  easily 

■  3  2 

with  the  aid  of  the  formulas  7.7,  7.22,  7.3;<  and  7.42. 


T.  37. 


j,  1  1  4  *  63 

*2(*> =  •  t*  ai s,n  t+  rr5 


cos  *  - 


3  19  ♦ 

-  (  — +  —  cos  *)  111  (  1  +  sin—)  + 


+  ( 4-  "  1tcos  *)  ln  sin  cos  *  L  ( *) 

5  2  5  *-0 


Ts  K  <*>  +  Js  E  (*> 


(8.0) 


The  derivatives  of  $  (♦)  can  now  be  obtained  from  the  formulas  8.6.  7.36, 
.7.45  -  7.48. 


t<  ^  2(  *)  g2  1 2  *  2 

- =  -  sin  *  -  - cos  —  -  — 

d*  125  25 


4r 

cot  f 


2  25  1+sinf 


+  (  -fj  sin  *  -  —  cot  y>  ln  (  1  +  5in  Y* 


34  3  ^  f 

+  (  —  sin  *  -  j  tg  j)  ln  sin  y 


-f-  si‘>  *  L  (*)  ♦  —  tg  -|"K  (*)  +  —  cot  *  E  ( *) 


(8.7) 


0"#  .,(*) 


8  *  202  ^  1  2  ♦ 

—  sin  —  +  -  cos  *  +  -  cosec  —  - 


d  t 


125 


25 


—  cosec  —  -  — 
<0  2  25 


1  ♦  sin  ~ 


T49  *  3  3  2  t  1 

[it  C0S  *  -  T  •  To  iec  T J 

f49  *  3  3  2  ♦  "1  ,  .  .  . 

[25  C0!  *  To"  co,ec  tJ1”  <’  +  ,inT) 


ln  sin 


-  }coi  +  l  (♦)  (sec2-!*  1  )  K  (*) 


^2  cosec2*  ♦  ij  E(  ♦) 


(8.8) 
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* 


Vic  still  have  to  evaluate  the  derivatives  of  the  angular  distance  V  with  respect 

to  \q.  These  derivatives  can  in  principle  be  obtained  from  the  spherical  cosine  theorem 

but  are  more  easily  found  from  three  fundamental  formulas  of  spherical  trigonometry. 

First,  we  make  the  change  of  variables  0  =  90°  -  6  ;  0  =  90°  -  0, 

0  0 

where  0  is  the  reduced  latitude  of  the  point  in  question.  The  three  fundamental  formulas 
are  (cf.  fig.  8. 1). 


sin 

* 

sin  A  - 
0 

cos 

8 

s  i  n  (  K  - 

V 

(8.  9) 

sin 

* 

c  o  s  A  = 
0 

sin 

0 

cos  0 

0 

cos 

0 

sin 

B  cos(\-X  ) 

0  0 

(8.10) 

cos 

♦ 

-  sin  3 

0 

sin 

0 

COS  0 

0 

cos 

0 

cos 

o 

* 

• 

(8.11) 

P 


Fig.  8.1 

P  *  the  north  pole  . 

=  the  fixed  point 

M  *  the  running  point 

A  «  the  azimuth 
o. 

The  angle  M^PM  is  equal  to  X-\q. 


Differentiating  (8.11)  with  respect  to  Kq  and  taking  (8.9)  into  account  we  obtain 


d  * 

- —  -  -  sin  A  cos  0 
ok  oo 

0 

Differentiating  once  more: 


8A0 


d\> 


(8.12) 


cos  0  cos  A 
o  o 


(8.13) 
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3A 

The  derivative  can  be  obtained  from  (8.9)  and  (8.10)  after  differentiation: 

8X 

o 

av  dAo 

cos  y  sin  A  —  +  cos  A  sin  y  -rr-  =  -  cos  0  cos  (X  -  A  ) 

o  a x  o  ax  o 

o  o 


a  y 


a  a 


cos  y  cos  A  — —  -  sin  A  sin  y  7—  =  -  cos  0  sin  8  sin(X-X  ) 

0  u  0  ax  o  o 

0  0 


From  this  we  have 
d  A 


sin 


a  a  f  “| 

y  — — —  =  cos  0  sin  fi  sin  A  sin  (X  -  X  )-  cos  A  cos  (X  -  A  ) 

ax  L  0  °  0  0  °j 


But  we  also  have: 


cos  (360  -A)  =  -  cos  A  cos  (X  -X  )  +  sin  A  cos  (X  -  X  )  sin  0 

OOO  00 


Thus 


aA  cos  a  cos  a 
_ 0  _  _ 

ax' 


sin  y 


When  we  substitute  this  expression  into'  (8.13),  we  obtain: 


a  y 


cos  a  cos  0  cos  A  cos  A 

_ o _ o _ 

sin  y 


(8.  M) 


Regarding  the  function  f  ^  as  a  function  of  the  four  variables  0Q,  0,  Aq,  A  we  obtain 

$„(*;  8  .  a  .  A  ,  A)  =  f,(y)  sin2  A  cos20  - 
3  o  0  1  00 


*  f„(  *)  cos  A  cos  A  cot  0  cos  0 
2  0  0 


where 


f,<*) 


f2<*> 


d  §2<*> 
d*2 
I 


sin  t 


dt 


(8.  15) 


(8.16) 


(8. 17) 


The  expretilons  for  and  aie  determined  further  by  the  formulas  (8.7)  and 
(8.8).  Thus  we  can  write  for  the  function  Z: 


Z  (*;  8  ,5.  A  ,  A)  =  S  ( *)  +  2  a  f(  *)  +  2q  #(*) 

0  0 

-  2a  [f  (♦)  sin^A  cos^S  -  f.(  ♦)  cos  A  cosAcosB  cos  fll  (8.18) 

1_2  oo2  o  °J 

where  S(+)  is  Stokes'  function  and  the  functions  !?(*).  ((♦),  f  (♦),  f  ( ♦)  are  given 

2  12 

by  (8.4),  (8.5),  (8.16).  (8.17).  The  quantity  ef  is  replaced  by  the  approximation  2  a. 

The  function  Z,  as  given  by  (8.18),  is  still  apparently  a  function  of  the  variables  A  and 
8  of  the  running  point.  In  reality,  Z  is  a  function  of  the  angles  +  and  A^  only  if 
we  assume  8q  to  be  given. 

From  the  spherical  triangle  rM^M  we  have: 

cos  8  cos  A  =  sin  *  sin  8  -  cos  8  cos  *  cos  A  (8.19) 

o  o  o 

Introducing  (8.19)  into  (8.18)  we  obtain 


Z(  «0:  *.  Ao)  =  S  (*)  +  2 a  f(*)  +  2q^(*)  - 


-  2a 


r  2  2 

f,(t)  sin  A  cos  8  -  f„(+)  cos  8  cos  A 

|_  1  o  o  2'  o 


(  sin  ♦  sin  8  *  cos  8  cos  t  cos  A  ) 

oo  o 


0 


(8.20) 


(o) 


Thus  we  can  write  the  expression  for  the  height  difference  N  at  the  fixed 
point,  determined  by  (8qi  Aq).  in  the  form: 


it  2n 

—  )  )  Ag*Z  (8  ;  A  )  sin  *  d  t 


d  A 


(8. 21) 


o 


Using  the  variables  *,  A  the  surface  element  do  =  sintdt  dA  . 

°  o  o 

The  quantity  A  g*'  is  connected  to  the  gravity  anomaly  Ag  by  the  equation 

2 

*■  i  /  2  2  6  1  2 

Ag  -  ul  *  ef  sin  {  Ag  ‘  Ag  (1  4  —  sin  8  ♦...)  * 


*  A g  ( 1  +  a  sin  8) 


(8.22) 
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where  sin  6  is  given  in  terms  of  the  coordinates  of  the  fixed  point  and  the  angles 

*  and  A  by  the  relation 
o  * 

2  2  2 

sin  6  =  sin  S  cos  +  +  2  sin  6  cos  0  sin  *  cos  *  cos  A 
o  oo  o 


The  integration  in  (8.21)  is  performed  on  the  circumscribed  sphere  with  the 
greacer  half-axis  as  radius.  Our  problem  is  now  solved. 


(8.23) 
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CONCLUSIONS. 


We  may  now  write  equation  (8.21)  in  the  form 


a  2  fl 


o  0 


N(0)  =  \  Ag*z(0)(  +  ;  8  n  >  A  ^  ,  8  .  A  dtdAJ 

2  n  y  o  o 


(9.  I) 


The  function  Z  is  given  by  the  expression 


(o)  r  22 

Zv  '(8  ,  ♦,  A  )  =  F(+)  +  2a  F  (*)  -  Fit)  sin  A  cos  8 
oo  L  o  1  o  o 


+  F_*  cos  8  cos  A  (sin  +  sin  8  -  cos  8  cos  *  cos  A  ) 

2  oo  o  o  o 


0 


+  2  q  Fs(*) 


There  we  have 


(9.2) 


F(  ♦)  =  —sin  *  S(t) 


F0(*)  ■  -j-  «»■•*  f(*> 


F,(+)  •  —sin  ♦  f,(+) 


Fa<  *>  =  —  sin  tfs(t) 


F3(*)  =  -5- sin  ♦  V'(  *) 

F(*)  is  sometimes  called  Helmert's  function.  Thus  we  have  as  our  final 
formula  for  N<0): 

11  2a 

N(°)  s  $  $  Ag  ["»  ♦  «sin2al  |f(  ♦)  ♦  2a  |f  ( 
2a  y  0  o  *-  -s,  L  V° 


(9.3) 

(9.4) 

(9.  8) 

(9.6) 

(9.T) 


2  2 

F^*)  sin  A^coi  8q  +  F^C  cos  A^cos8o  (sin  *  sin  6^- 


-  cos  A  cos  8  cos 
0  o 


♦)|+  2q  Fg(  ♦  )J 


d  td  A 


(9.8) 
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or  with  the  same  degree  of  approximation: 


fl  2  fl 


where 


N(o)  =  —  \  $  Ag  2<0)<flo:  *'  V  d*dA0 

■  2 » y  o  o 

2(0)(8o;  Aq)  =  F(*)  £l  +  a  sin2ej  +  2a  j^M*)  - 

2  2 

-  F . ( ♦)  cos  fl  sin  A  + 

1  o  o 

1  2  2 
+  —  F^(+)  s*n  ®®0s^n  *  cos  A0  *  F2(*)  COS  do  C°S  *cos  A0 


(9.  9> 


>] 


♦  2q  F3(t) 

2 

Here  sin  0  is  given  by  (8.23)  as  a  function  of  0  q,  *.  A  .  The  quantities  in 

(9.9)  and  (9.10)  are  defined  as  follcws: 


(9.10) 


=  the  height  difference  between  the  geoid  and  the  ellipsoid 
a  -  the  greau:  half- axis  of  the  ellipsoid 

q  =  the  quotient  between  centrifugal  force  and  gravitational  force  at  the  equator 

y  -  the  gravity  of  the  theoretical  earth 

Bo  =  reduced  latitude  of  the  fixed  point 

=  angular  distance  between  fixed  point  and  running  points. 

The  anomaly  Ag  in  (9.9)  is  regarded  as  a  function  of  the  reduced  latitude  and  the 

longitude.  The  functions  F,  F ^ ,  F  ^  .  F^,  F^  occurring  'n  t*ie  expression  for 

(0  ;  +  ,  A  )  are  given  by  the  formulas  (  cf.  9.3  -  9.7,  8.4,  8.5,  8. 16  and 
o  o 

8. 17). 


F(  ♦) 


sin 


coscc 


3 


3  * 

—  cos  t  In  (sin  —  ♦ 


5 

Y  cos 


* 


(9.  11) 
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F  (+)  =  sin  <t 
o 


.1 1 


Ft-t 


*  42  *  164 

C0S6C  T  +  25  Slr  ~+  TIT  C0S  ** 


1  1 


<50  cos  *  ■  fo>  ln  +  sin  T>  +  (  —  cos  *  +  — )  in  sin-r- 


50 


1  0' 


4  To  cos  *  E(+)] 


Fi(t)  =  sin  *  [(  sin  ^  C0S  *  “  To  •  — ] + 

l+sinT 


,  49  3  3  1  *  o 

4  (  To-  cos  *  4  To  4  2°'“7+,‘  Ln  (1+  sin  ~T)  '  To‘°**U*) 

S1"  2 


*2l  E(t)]  ‘IT  cos  7  4  2l[sin  TK(t)(cos  T 

C0t4  -  2  •  E(  +) — 1+(—  COS  — 

2  v  ’  sin  *  J  1  25  2 


* 

c°sT 


)  ♦ 


sin  —In  sin 


3  *  3  1 

cos  t  -  —  cos  —  -  —  —  .  ) 

5  210  *  ’ 

cosT 


F2(#)  =  ih  ‘ln  *  '  h  C0S  T  ’  Jl - cot  T"  4 

Usin-j 

4  (7TsIn  *  •  TT  cot  T)  ln  (  1  4  sin  7*  4  *in  *• 

*  To  tgT)  ln  sin  T '  To  $in  *  1  (*>  +  Is  T"k^ +iTcot* 


f3<*)  =Ysint(-1  +  2  sin  4+  co»  *+  (3-5  cos  *)  In  sin  -1- 


♦  ■  in-s-)  +  3  cos  *  L(  ♦) 


(9.12) 


(9.13) 


E(t) 

(9.  M) 


-(3  +  5  cos  *)  in  ( 1 


(9.15) 


t 
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In  che  last  four  formulas  K(+)  and  E(+)  are  complete  elliptic  integrals  of  the 


first  and  second  kind  and  of  modulus  cos  —  . 

Thus  , 


a 

"2 


K(*) 


-  1 


d  <f 


O  dr  O 

1  -  cos  —  sin  if 


(9.16) 


E<*)  =  \ 


2  |  2 

I  -  cos  —  sin  if  d  if 


(9.17) 


L  (  +  )  is  a  new  function,  given  by  the  formula: 


L(*)  =  S  ( 


1 


0  ^i~~ 


-  — )  dx 
X 


2x  cos  *  +  x 


(9. 18) 


(cf.  formulas  7.34  et  scq.) 
This  function  can  be  tabulated  (see  table  below)  for  some  values  by  the  use  of  its 
derivative 


d  L(  *)  *  ,  *  *  ,  *. 

-d*”“  tg  Tln  sinT  *  cotTla  <’  ♦  slnT) 


(9.19) 


and  the  fact  that  L  (o)  *  -g-  . 

After  evaluating  according  to  formula  (9. 9)  (in  principle  for  every  point  of  the 

earth)  we  must  compute  the  correction  This  can  be  done  by  replacing  A  g* 

2  2N^a^v  2 

(=  Ag  +  Ag  sin  8)  in  (9.8)  or  (9.9)  by  -  ■■  ^ ‘  cos  0  (cf.  eq.  3.9  * 

3.11  and  1.1). 

It  will  be  sufficient  to  take  into  consideration  just  the  term  F(t)  from  (9.8)  or 
(9.9).  As  we  neglect  terms  of  order  we  do  not  need  more  than  this  term  of  our  expression 
for  0N(,). 


Thus, 


where 


a  2a 


\  S  Ag(1)  F(t)  d*dA 


2a  y  o  o 


(1)  2N(0)y  2. 

=  - L  cos  0 


(9. 20) 


(9.21) 
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t 


We  have  now  solved  our  main  problem,  and  it  remains  only  to  test  the  results 
on  a  suitable  model.  The  solution  is  given  by  the  formulas  (9.9)  and  (9.20). 


N 


+  a  N 


(1) 


(9. 22) 


The  new  functions  are  compared  to  Sagrebin's  functions  in  a  diagram  below. 

A  special  test  of  the  function  <p  has  been  performed,  by  directly  summing  the  series 
on  an  electronic  computer 


2 


OO 

=  £ 
n=  2 


1 

(n-1)2 


P 

n 


(cos 


+) 


for  different  values  of  +.  The  results  obtained  agree  completely  with  those  found  from 
the  explicit  formula  (7.33)  for  t>2(+). 
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THE  TEST  MODEL. 


For  a  simple  test  of  the  resolvent  function  Z,  the  following  "model  earth"  is 
designed  (Fig.  10.1), 

The  theoretical  gravity  y  is  obtained  horn  the  requisite  massdistribution  to  make  the 
ellipsoid  an  equipotential  surface.  The  "formula  of  Bruns"  then  holds  for  y: 


1  dy  1  (  1 

r  an  (pm  pn 


2  w 


(10.1) 


To  obtain  a  disturbance  potential  T  with  the  gravity  centres  of  the  theoretical 

earth  and  the  actual  earth  coinciding  we  put  four  unit  masses  into  the  ellipsoid  in  the 
o 

plane  X  =  0,  180  .  Furthermore,  we  put  a  negative  mass  of  -4  units  in  the  centre. 
The  radius  of  the  circumscribed  sphere  and  the  gravitational  constant  are  taken  to  be 
the  unity. 


Thus 


4 

T  =  E 
i=1 


(10.2) 


oj 


_L  ±JL 

y  d  n 


T 


(10.3) 


Vve  know  the  disturbance  potential  a  priori,  and  therefore  we  can  simplify  the 
computations  to  one  single  step.  This  means  we  have  to  use  the  boundary  condition 


where 


(1  + 


.  2T 
q)_  + 


8  T 

8p 


T  =  the  total  disturbance  potential 
Ag"*=  the  gravity  anomaly  for  our  resolvent. 

Furthermore  we  have  the  relation 


Ag  =  A  g~(  1  +  2a  sin  6) 


(10.4) 


(10.5) 
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Substituting  this  A  g  into  the  resolvent  formula  we  obtain 


a  2a 

T  =  2“  i  $  AgZ  d ♦ d A 
o  o 


(10.6) 


This  value  should  coincide  with  the  value  from  the  Newtonian  computation. 
The  computation  of  A  g. 

Starting  from  (10.2)  we  have: 


Symbols 


,  T  4  4  (nr..)  (nr  ) 

il.  E±±.  4f  £  — lL+4_il  (10.1) 

dn  i=1  dn  rij  dn  roj  i=1  r  3  r3  . 

ij  oj 


oj 

roi 


the  vector  from  the  origin  to  a  point  on  the  ellipsoid 
the  vector  from  the  origin  to  a  disturbing  mass 
die  length  of  a  vector 


'ij  *  *  r' 


io  jo 

s  the  outer  normal  of  the  ellipsoid  (unit  vector). 

*  unit  vector  ^ 

Thus  we  have  to  investigate  the  quantities  (u  -  r^)  and  (n  •  ~Qj)- 

The  normal  of  the  ellipsoid  is  given  by  the  formula 


n 

r 


"  •  “oj  +  “  2  »j  rj 


(10.8) 


where  Bj  is  the  unitvector  of  the  tangent  (towards  the  north  pole)  of  the  sphere  at  the 
reduced  point  corresponding  to  j. 

Thus 


(n'roj)  ■  r0j  +  *W  »0j>  "  roj+0(a  > 


as  (£}-  7  )  ■  0(«) 


(10.8) 


(n 


•7oi)  *  <V  -7oi  “<W  + 


•,ol,1“Mj(*j,,o|) 
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(r  .■  t  .)  =  r  .  cos  v,  . 
oj  01  01  i] 


(10.11) 


where  v,,  is  the  angle  between  the  vectors  7  .  and  7  .. 
ij  B  oi  oj 

The  quantity  (^  •  7q  , )  remains  to  be  found.  In  Cartesian  coordinates 
(  x  ,  y ,  z)  we  have 


Here 


r  =  TT  cos  6,  +  z  sin  6.  +  0(a) 
oi  i  i  i 


B.  =  -R.  sin  fl  .  +  z  cos  B  .  +  O(a) 
J  J  J  J 


R  «  x  cos  X  +  y  sin  X 


(10.12) 

(10.13) 

(10.14) 


(8.  r  .)  =  -cos  g.  sin  8  .cos  (X.-X  )  +  sinfl.  cosfl.  +  0(a) 

J  oi  i  j  J  *  i  ) 

We  do  not  want  the  quantities  fl^  and  V  to  enter  into  the  final  result,  since  we 
integrate  over  *  and  A  (*A^).  Here  B^.  X^  and  6^  X^  are  given  beforehand 
and  may  enter  explicitly  into  the  formulas  for  Ag  and  T. 

We  have  the  formulas 


sin  .  sin  A  *  cos  0.(sinX.cosX,  -  sin  X,  cosX.)  (10.15) 

jk  jk  J  j  k  k  j  . 


cos  ♦  =  sin  6  .  sin  8.  ♦  cos  fl  cos  fl.  (cos  X.  cos  X  ♦ 

)k  jk  J  k  jk 

sin  X,  sin  X,  )• 
j  k 

(cf.  eq.  (0.9)  and  (0.11)). 

From  the  equations  (10.15)  and  (10.16)  we  can  solve  for  the  quantities 


(10.16) 


cos  0 j  cos  X j  (■  Cj)  and  cos  0^  sin  X^  (*s^). 


cos  t  -  sin  fl,  sin  fl  . 
ik  k  1 

cos  fl. 


cos  X^  -  sin  t  sin  A^sin  Xfc  (10.  11) 
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cos  V  -  sin  0  sin  0 

s.=  - - - - — - -  sin  X,  +  sin  sin  A  .  cos  X  (10.18) 

j  cos  8^  k  jk  jk  k 

A  j  ^  =  A  and  *jk  =  *  ve  integration  variables. 

Furthermore 


cos  + 


sin  0 ,  sin  0  . 
k _ 1 

COS  0  , 
k 


sin  *  cos  A  -  sin  0  ^  cos  0^ 


sin  0, 


(cf.  (8.  10)  and  (8.  11)) 

This  expression  can  be  used  when  cos  0^  approaches  zero. 

By  permutating  the  angles 

90°  -  0  .-*•  * -*-90°  -  0 
)  k 

X  -  X^— »A  in  (10. 16)  we  obtain 

sin  0  .  »  sin  0  ^  c0'  *  +  cos  0 ^  sin  #  cos  A  (10.19) 


Since 


-90°  <  0  .  <  90° 
J 


we  have  .simply 


The ')  educed” 


cos  0 j  =  yi  -  sin2©^  (always  positive) 

angle  between  the  vectors  7  ,  and  T  .  is  given  by 
lj  o  i  o  j  °  ' 


(10.20) 


cos  *  sin  04  sin  0 ^  +  cos  0J  (c^  cos  X(  ♦  8^  sin  Xj)  (10.21) 

To  obtain  the  physical  angle  v^  in  (10. 11).  we  have  to  replace  0^  by 
a 

(  flj  •  —  sin  20 1)  and  analogously  for  0^. 

Thus  by  series  expansion 


cos  v 


U 


2  2 

cos  +ij(1+  a(sin  0^+sin  0.))  -  BasinB^sinB^ 


(10.22) 
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Expression  (10.15)  reduces  to 


cos  8  .stn  8  . 

(I .  •  7  .)  =  -  - -  J  (c.cos  X+s.sin  X.)  +  sin  8, cos  8, 

J  oi  cosdj  j  i  ]  i  i  j 

Thus 

-  —  2 
a  sin  2  8  .(8  .  •  r  .)  =  2asin8.(sln8,cos  8.-  cos 6, sin  8. 

J  J  oi  J  i  J  i  j 


( c  cos  X.  +  s  .sin  X.)) 
j  i  J  i 


(1  0. 23) 


We  have  also 


2  2  9 

r  =  r  .  +  r  ,  -  2t  .  t  cos  v . 
y  oi  oj  oi  oj  ij 


(10. 24) 


Prom  the  equation  of  the  ellipsoid  we  have 


2  2 
r  .  =  1  -  2a  sin  8. 
oj  j 


(10.25) 


2  2  z 

X  +  y  + — n  =  i 

1  -  e 

( z  axis  along  the  axis  of  rotation) 

where  e  is  the  first  eccentricity  of  the  ellipsoid. 


2  2  2  ,  „  2  „  2v 
x  +y  ♦  z  =  1  -  2a  z  +  O (a  ) 


Now  we  can  write  (cf.  10.7) 


z=  sin  8  j  +  0(  a) 


il.  -  lu  .4^21 

On  i=1  r  ®  r  3 

ij  oj 


4  4 

)  *  -  -  E  - (r  ,-r  .  cos  v ,  - 

,  ,«  I-  ,,  ,3  °J  01  “ 

oj  ij 


2 

2a  sin  8. (sin  8,  cos  8  -  cos  8,sin  8,  d..)) 

j  i  J  i  j  j  i  ' 


( 1  0. 26) 


where 


d^  •  Cj  co,  Xj  +  ,j  sin  Xt 


The  formula,  given  here  contain  all  the  information  necessary  to  test  the  modified  t 
Sagrebin  resolvent  by  the  aid  of  an  electronic  computer. 


le  position  of  the  five  pointmasscs  in  the  test  model.  The  central  mass  is  equal 
-4  and  the  remaining  masses  are  equal  to  +1 . 
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RESULTS  FROM  THE  TEST  MODEL. 


Test 

N  :  o  1 

a  -  0 

0 

0  ^  60 

X  =  1  5° 

Compartment: 

da  =  15° 

0 

dX  =  1  5 

Newtonian  T  = 

0.81870 

Resolvent  T  = 

0. 84866 

( Stokes) 

Error 

+0.02996 

Test 

N  :  o  2 

a  -  0 

8  =  60° 

11 

cn 

O 

Compartment: 

'da  =  5° 

d\  =  5° 

Newtonian  T  = 

0.81870 

Resolvent  T  = 

0.81790 

( Stokes  ) 

Error 

-0.00080 

Test  N :  o  3 


a  =  0.001 
Compartment: 

Newtonian  T  = 
Resolvent  T  = 
Error 


o 

e  =  60 

o 

da  =  15 

0. 81962 
0,85102 
+0. 03140 


X  = 
dX  = 


o 


Test  N  :  o  4 


o  =  0.001 

Compartment: 

Newtonian  T  = 
Resolvent  T  = 


8  =  60° 
da  =  5° 

0. 81962 
0. 81998 


X=1  5 


o 


d 


_  o 

D 


Error 


f 0. 00036 
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Couclusions  from  numerical  tests. 


It  has  been  possible  to  verify  tils'  principal  formula  v2(x)  by  the  aid  of  a 

numerical  series  expansion  on  an  electronic  computer.  This  calculation  has  proofed 

that  the  new  derivation  is  correct.  The  final  resolvents  includes  a  number  of  functions 

which  all  are  derived  from  this  principal  formula.  The  whole  assembly  of  formulas  is 

verified  by  the  aid  of  the  test  model.  The  study  of  this  model  has  shown  that  error  in 

the  computations  is  decreasing  with  approximately  forth  power  of  size  of  the  integra- 

o  o 

tion  compartments.  .Already  with  a  compartment  01  1x1  the  resolvent  solution  will 

fit  better  than  10"®  of  T.  The  new  resolvent  tor  the  ellipsoid  has  given  an  error  which 
was  only  half  of  that  corresponding  to  the  same  solution  according  to  Stokes  for  a  real 
sphere  (  5°  x  5°  ).  One  can  expect  both  formulas  to  have  approximately  the  same 


accuracy. 


COMPARISON  BETWEEN  SAG  REBINS  FUNCTIONS  AND  THE 


NEW  FUNCTIONS. 


SagreJjin: 


4  1  1/  12  *  253 

F  (  ♦)  =  sin  *  (  —  -  —  cosc-c  —  +  —  sin  —  +  — —  cos  *  + 

o  I  o  1  2  25  2  250 

3  *  2  *  11  *  2  * 

+  —  In  (sin— +  sin  —  )  +  —  cos*  X  n  ( s  i  n  — +  sin  —  ) 


+  To*  <*>  ’Is  E  <*> 


Corrected  F  (V)  function: 
o 


Fo(*)  =  sin  *  (  -y-  -1  cosec  ii  sinA  + 


164  11  3  * 

+  -TTrr  cos  +  •  (—  cos  *  -  —  >  ln  (’  +  sin  — )  + 


1  25 


1  0 


♦  <73  cos  *  *  To*  lD  sin  T*  To  cos  *  L(t)  +  JoK(*)'IsE(*) 


In  a  similar  way  the  three  remaining  functions  F.(+)  F  (*)  and  F  (+)  are 

»  &  « 

changed. 


Sagrebin: 


«'2(  +  )  =  J-  -  2  *in~+  2  sin2  yin  (‘in  y  +  sin2  ^~) 


Correct  function: 


*  * 

?2<  ♦)■-!♦  cos  *  ♦  2  sin  y  ♦  (1-  cos  *)  In  sin  —  - 


(  1  ♦  cos  *)  In  (  1  ♦  siny)  -  cos  ♦  J<2^ 


where 


$  4 
0 


Y  (  1  •  w*  ♦  r(x)) 


dx 


-S  crk-4-; 


r(x)  x 


•)  in  x  dx 


etc . 


From  a  theoretical  point  of  view  the  new  results  seem  to  be  of  the  greatest 
importance  since  hitherto  modern  geodesy  is  based  on  the  fact  that  the  errors  in  the 
generally  used  Stokes  approach  could  be  determined  according  to  Pizetti  and  then 
finally  according  to  Sagrebin.  However,  Sagrebin  has  shown  that  Pizetti's  method  is 
insufficient  to  give  any  real  information,  Sagrebin's  own  theoretical  study  is  probably 
one  of  the  most  important  in  modern  geodesy,  in  spite  of  the  fact  that  it  is  incorrect 
in  all  new  functions. 

It  is  supposed  that  some  of  the  fundamental  -questions  in  geodesy  can  be  solved  by 
the  aid  of  this  new  study. 


0 


TABLE 

OF  FIJNC 

T  I  0  N  S 

(15°) 

X 

F(X) 

FO(X) 

F1(X) 

F2(X) 

F3(X) 

0 

0 

1.0000 

-.5000 

.0000 

-.5000 

.0000 

15 

1.1121 

-.3618 

-.1031 

-.7901 

.8928 

30 

.4736 

-.0683 

.1625 

-.7924 

.7923 

45 

-.3070 

.1793 

.5495 

-.6429 

.1359 

60 

-.8957 

.3103 

.8677 

-.4053 

-.6297 

75 

-1.1049 

.3137 

.9841 

-.1365 

-1.1127 

90 

-.9142 

.2197 

.8598 

.1128 

-1.1150 

105 

-.4500 

.0808 

.5531 

.3025 

-.6825 

120 

.0773 

-.0487 

.1872 

.4066 

-.0528 

135 

.4578 

-.1281 

-.1034 

.4152 

.4705 

150 

.5590 

-.1391 

-.2275 

.3347 

.6590 

165 

.3702 

-.0879 

-.1713 

.1859 

.4574 

180 

.0000 

.0000 

.0000 

.0000 

.0000 

F(  *)  = 

Y s*n  *  S  (  +) 

i  i 

=  Tsin*  47 

CO 

7  2n+ 1 

E  - -  P 

_  n*  1 

n=2 

„(c°s  ♦) 

F0(+)  =  Tsin  *  =  T,in  *  <*•,(♦)  +  2  e>2  (*)) 


T  sln  y 


oo  2 

g  n  (  2n+  1 ) 
n  =  2  (n-  1  )2(2n+3) 


P  (cos  ♦) 

o 


CO 

♦  2  E 
n=2 


(2n-f1) 

(n-  1)2(2i»3) 


P  (COS 
n 


F,(+)  =  ~  sin  *  f,(*)  =  ~  sin  * 


d  «>2(*) 

“77“ 


f2(+) 


Y  sin  *  f2(  *) 


1  d  0  ( t) 

2  dt 


F 


3<*> 


2 


sin  t  $T(  *)  =  -1  sin 


GO 

♦  E 
n  =  2 


(2n+1) 

(n-1)2 


PQ  (cos  *). 


TABLE 

0  F  F  U 

X 

F(X) 

F0(X) 

0 

0 

1.0000 

-.5000 

5 

1.2165 

-.5055 

10 

1.2146 

-.4456 

15 

1.1121 

-.3618 

20 

.9410 

-.2663 

25 

.7224 

-.1666 

30 

.4736 

-.0683 

35 

.2100 

.0242 

40 

-.0543 

.1076 

45 

-.3070 

.1793 

50 

-.5371 

.2375 

55 

-.7357 

.2813 

60 

-.8957 

.3103 

65 

-1.0121 

.3248 

70 

-1.0820 

.3255 

75 

-1.1049 

.3137 

80 

-1.0821 

.2908 

85 

-1.0169 

.2588 

90 

-.9142 

.2197 

95 

-.7804 

.1756 

100 

-.6230 

.1286 

105 

-.4500 

.0808 

110 

-.2700 

.0342 

115 

-.0916 

-.0095 

120 

.0773 

-.0487 

125 

.2292 

-.0822 

130 

.3577 

-.1068 

135 

.4578 

-.1281 

140 

.5258 

-.1396 

145 

.5596 

-.1432 

150 

.5590 

-.1391 

155 

.5251 

-.1281 

160 

.4608 

-.1106 

165 

.3702 

-.0879 

170 

.2589 

-.0610 

175 

.1331 

-.0312 

180 

.0000 

.0000 

T  I  0  N  S 

<«*> 

F1(X) 

F2(X) 

F3(X) 

.0000 

-.5000 

.0000 

-.1071 

-.6567 

.4613 

-.1288 

-.7409 

.7364 

-.1031 

-.7901 

.8928 

-.0407 

-.8123 

.9460 

.0504 

-.8120 

.9080 

.1625 

-.7924 

.7923 

.2881 

-.7561 

.6140 

.4195 

-.7055 

.3894 

.5495 

-.6429 

.1359 

.6715 

-.5706 

-.1292 

.7792 

-.4907 

-.3895 

.8677 

-.4053 

-.6297 

.9330 

-.3165 

-.8369 

.9723 

-.2262 

-1.0004 

.9841 

-.1365 

-1.1127 

.9684 

-.0490 

-1.1696 

.9262 

.0346 

-1.1696 

.8598 

.1128 

-1.1150 

.7726 

.1843 

-1.0103 

.6688 

.2478 

-.8630 

.5531 

.3025 

-.6825 

.4307 

.3475 

-.4794 

.3070 

.3823 

-.2656 

.1872 

.4066 

-.0528 

.0761 

.4201 

.1475 

-.0219 

.4228 

.3249 

-.1034 

.4152 

.4705 

-.1657 

.3975 

.5773 

-.2073 

.3704 

.6409 

-.2275 

.3347 

.6590 

-.2271 

.2913 

.6320 

-.2075 

.2413 

.5631 

-.1713 

.1859 

.4574 

-.1219 

.1263 

.3222 

-.0633 

.0638 

.1664 

.0000 

.0000 

.0000 

732501  -  514  A4  -  Ixlrmn 


732501  -  514  A4-  1x1mm 


73  2501  -  514A4-  1*1 
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APPENDIX. 


TABLE  FOR  THE  RESOLVENT  OF  AN  ELLIPSOIDAL 
SURFACE  OF  REFERENCE. 


The  table  can  be  used  for  determining  the  height  over  the  geoid  in  the  following 
manner. 

We  start  with  formula 


t  2lt 


-  ■  2  2 

=  -r— —  ^  $  A  g  1  +  (sin  6  cos  *  +  2sin8cos8  sin*  cos* 

2  s  y  2  o  oo 

oo  *“ 


cos  A  +  cos  ^6  sin2*  cos2A)  (F(*)  +  i2F  (  *)-i2  F  (  *)  s in2A  cos28  + 
O  J  0  1  o 

2  2  2  2 
+  i  F  (  *)  cos  A  cos  6  sin  *  sin  8  -  i  F  ( *)  cos  A  cos  8  cos  *  + 

Z  0  0  Z  O 


+  2q  F g(  ♦))  d  *d A 


This  formula  is  transcribed 


H 


where 


3  p  -» 

N(0)  -  C  $  \  Ag  ♦  corrj  (E  a^f*))  dS  [l3.lj 


corr^  =  O  (  i  ) 


F.,  (♦)  =  F  (*)  def. 


a 

v 


v  =  0,  1 ,  2,  3 
v  =  -1 


Every  a^ 


is  a  function  of 


8  =  reduced  latitude  for  the  fixed  point, 

o 

*  «=  angle  between  the  fixed  and  the  running  point 

A  =  azimuth  for  the  running  point  (measured  in  the  onhonormal 

system  determined  by  the  fixed  point). 


We  substitute  now  for  equivalent: 

N(0)  =  C  i  S  ig  F(*)  [~l^f  ' 

where 


+  corr^  dS 


F  (  ♦)  =  F  (  *)  is  independant  of  the  position  of  the  fixed  point.  The 

‘  ’  2 
factor  (1  +  corr^  E  ( ay  Fy) :  F(+)  =  1  +  G(0q;  *,  A)  •  O  (i  )  is  near  to  1 

for  any  point  on  the  ellipsoid  and  equals  1  for  every  point  on  a  spherical  surface.  This 

(O) 

fact  enables  us  to  use  the  following  approximation  formula  for  N 


N(0,=  CEE  Ag^  F(*l)-K(floj  +l.  A..)  A  t  j  [13.3J 

1  j 

where  the  magnitude  of  K  in  accordance  with  our  considerations  above  is  near  to  1. 

Now  an  appropriate  network  for  t.  and  A_  is  defined  and  we  have  in  our 
example  chosen: 


Aij  "  i  *  10 


i  =  0.1 . 36 

j  =  0 ,  1 . 3  6  for  every  i 


As  reduced  latitude  for  the  fixed  point  we  have  taken  0^  =  +  30  .  Now  we  have 

computed  the  values  of  K(30;  A  i .)  headed  under  +.  ,  furthermore  the  reduced 

latitude,  0  ,  for  (  V  ,  A  ^. )  in  the  orthonormal  system  determined  by  the  axis  of  the 

ellipsoid  and  the  difference  in  azimuth,  1  -  1  ,  between  the  fixed  and  the  running  points  in 

ij  o 

the  same  system. 

The  system  of  equations  for  the  coordinate  transformation  is: 


cos  0,,  cos  (1,,  -  1  j-  cos  #,  cos  0  +  sin  t  cos  A.. sin  0 

ij  ij  o  to  1  lj  o 

cos  0tj  sin  { 1 t j  -  1^=  sin  sin  A 

sin  0,  .  =  cos  *,  sin  6  -  sin  *.  cos  A, .  cos  0 

ij  i  0  i  ij  o 

The  table  is  constructed  so  that  for  any  point  where  Ag  is  given  one  can  get 
an  approximate  value  of  the  K- function  .  We  have  then  to  multiply  K  with  its  corre¬ 
sponding  value  of  R[t)  and  with  that  element  of  surface,  that  the  actual  point  is  in¬ 
tended  to  represent,  to  get  the  actual  partial  sum  of  |l3.3j  •  (Note  that  F  (  *)  ■  -y 
sin  t  S{  *),  where  S  (♦)  is  Stokes'  function,  and  that  therefore  "the  element  of 
surface''  lacks  the  sin  factor). 


# 

* 

Exampl 

1 
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e  of  tabulated  resolvent  values  according  to  (13.3). 

BETAs  Reduced  latitude. 

L-L(O)  = 

Longitude  difference. 

PS!-  5.0  BHTA 

L-L(O) 

Z(0)/F(c)  -  K(«a>  *t.  A,.) 

25. 0G  ° 

.00 

0  1.0035 

25.07 

.96 

1.0035 

25.2  9 

1.89 

1.0032 

25.64 

2.77 

1.0029 

26.12 

3.58 

1.0025 

P-6.72 

4 . 29 

1.0021 

27.41 

4.88 

1.0017 

28.18 

5.33 

1.0014 

29.01 

5.63 

1.0012 

29.87 

5.77 

1.0011 

30.74 

5.73 

1 .0012 

31.  BO 

5.52 

1.0013 

32.40 

5.13 

1 .0016 

33.14 

4.57 

1.0020 

33.77 

3.86 

1.0024 

34.30 

3.02 

1.0027 

34.68 

2.08 

1.0030 

34.92 

1.06 

1.0032 

3?.  00 

360.00 

1.0033 

, 

34.92 

358.94 

1 .0032 

34.68 

357.92 

1.0030 

34.30 

356.98 

1,0027 

33.77 

356.14 

1.0024 

33.14 

355.43 

1,0020 

32.40 

354.87 

1,0016 

31.60 

354.48 

1.0013 

30.74 

354.27 

1,0012 

29.87 

354.23 

1.0011 

29.01 

354.37 

1.0012 

28.18 

354.57 

1.0014 

27.41 

355.12 

1.0017 

26.72 

355.71 

1.0021 

26.1? 

356.42 

1.0025 

25.64 

357.23 

1.0029 

25.29 

358.11 

1.0032 

25.07 

359.04 

1.0035 

25.00 

.00 

1.0035 
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PSI- 


EFT  A 

L-L(O) 

Z(0)/F(0) 

20. 00 

.00 

1.0059 

20.14 

1.84 

1.0058 

20.55 

3.64 

1.0056 

21.23 

5.-34 

1.0052 

22.13 

6.92 

1.0047 

2.3.31 

8.33 

1.0043 

2b ,  06 

9.52 

1.0038 

26.17 

10.48 

1.0035 

27.79 

11.15 

1.0032 

29.50 

11.51 

1.0031 

31.23 

11.54 

1.0031 

32.95 

11.21 

1.0033 

3/l .  58 

10.52 

1.0036 

36.09 

9.47 

1 .0040 

37.42 

8.08 

1.0045 

38.51 

6.37 

1.0049 

39.32 

^.40 

1.0052 

39.83 

2.25 

1 .0054 

40.00 

360.00 

1.0055 

39.83 

357.75 

1.0054 

39.32 

355.60 

1.0052 

38.51 

353.63 

1.0049 

37.42 

351.92 

1.0045 

36.09 

350.53 

1 .0040 

34.58 

349.48 

1.0036 

32.95 

348.79 

1.0033 

31.23 

348.46 

1.0031 

29.50 

348.49 

1.0031 

27.79 

348.85 

1.0032 

26.17 

349.52 

1.0035 

24.66 

350.48 

1 .0038 

23.31 

351.67 

1.0043 

22.16 

353.08 

1.0047 

21.23 

354.66 

1.0052 

20.55 

356.36 

1.0056 

20.14 

358.16 

1.0058 

20.00 

.00 

1 .0059 
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BETA 

L-L(O) 

15.00 

o 

o 

• 

15.20 

2.67 

15.80 

5.28 

16.79 

7.77 

18.14 

10.08 

19.81 

12.17 

21.77 

13.97 

23.97 

15.44 

26.36 

16.53 

28.88 

17.19 

31.46 

17.39 

34.03 

17.07 

36.52 

16.19 

38.83 

14.75 

40.89 

12.71 

42.62 

10.13 

43.91 

7.06 

4^.72 

3.63 

45.00 

360.00 

44.72 

356.37 

43.91 

352.94 

42.62 

349.87 

40.89 

347.29 

38.83 

345.25 

36.52 

343.81 

34.03 

342.93 

31.46 

342.61 

28.88 

342.81 

26.36 

343.47 

23.97 

344.56 

21.77 

346.03 

19.81 

347.83 

18.14 

349.92 

16.79 

352.23 

15.80 

354.72 

15.20 

357.33 

15.00 

.00 

Z(0)/F(0) 

1 , 0080 
1 .0079 
1.0076 
1.0072 
1.0066 
1.0061 
1.0055 
1.0051 
1 .0048 
1.0046 
1.0047 
1.0049 
1.0053 
1.0058 
1.0063 
1.0068 
1 .0073 
1.0075 
1.0076 
1.0075 
1.0073 
1.0068 
1.0063 
1.0058 
1.0053 
1.0049 
1.0047 
1 .0046 
1 .0048 
1 .0051 
1.0055 
1.0061 
1.0066 
1.0072 
1.0076 
1.0079 
1.0080 


GRAVITATION 


Historical 

The  study  of  gravitation  has  attracted  an  exceptional  amount  of  interest  in  recent 
years.  The  foremost  reason  for  this  is  found  in  the  rapid  development  which  is  occuring 
at  present  in  the  disciplines  connected  with  geophysics.  The  problems  connected  with 
satellites  and  robot  weapons  have  accented  this  interest  further.  The  general  manner  of 
development  is  briefly  as  follows. 

In  Sweden  and  the  other  "western  countries”  developments  have  previously  followed 
the  classical  representation  of  STOKES.  In  the  application  of  this  theory  one  has  sought  to 
obtain  gravity  material  by  direct  measurement  of  gravity  in  various  places  over  the  whole 
earth  which  then  was  applied  to  determine  the  so  called  geoid.  The  auxilliary  surface 
can  be  said  to  constitute  the  gravitations  model  earth.  If  the  continents  were  rut  through 
by  a  network  of  canals,  the  mean  water  surface  thus  obtained  would  define  the  geoid.  One 
difficulty  with  the  application  of  STOKES’  formula  for  the  determination  of  the  geoid  is 
that  no  masses  can  lie  outside  the  actual  geoid.  In  the  practical  application  of  STOKES’ 
theories  one  is  therefore  forced  to  eliminate  by  some  artifice  the  masses  which  lie  outside 
die  geoid.  Western  scientists  have  not  been  successful  in  solving  this  problem.  At  the 
same  time,  many  of  the  works  which  have  been  carried  out  toward  this  goal  will  certainly 
be  significant  in  various  other  connections  even  if  they  do  not  succeed  in  the  solution  of  the 
main  problem.  To  correctly  understand  the  older  manner  of  consideration  it  can  be  useful, 
to  give  the  historical  development.  In  the  triangle  measurements  in. India,  J.  H.  PRATT  in 
t8S5  found  Urge  triangle  misclosures  which  were  measured  on  the  surface  of  the  earth  in 
the  vicinity  of  the  Himalaya  mountains.  From  this  tie  drew  the  conclusion  that  the  large 
masses  of  mountains  had  deep  "roots”  with  less  density  than  the  surrounding  area  so  that 
the  visible  outer  masses  were  completely  compensated  at  a  lower  level.  In  such  a  manner 
he  could  assume  that  two  equally  large  land  areas  always  contain  masses  of  equal  size. .  The 
assumption  for  this  is  only  that  the  two  volumes  compared  include  all  the  masses  down  to 
the  earths  inner  surface  of  compensation.  The  hypotheses  presented  by  PRATT  have  led  to 
a  special  uiscipline  which  is  called  "isostasy".  Of  great  significance  for  the  application  of 
isostatic  hypotheses  wu  the  determination  of  the  "depth  of  compensation”  which  is  required 
to  obtain  equilibrium  among  the  masses.  On  this  point  views  differ.  In  his  own  theories 


PRATT  assumed  that  the  depth  of  compensation  was  everywhere  constant  and  many  writers 
said  that  this  depth  should  be  measured  from  the  unknown  surface  of  the  geoid.  The  Ame¬ 
rican  JOHN  F.  HAYFORD  was  the  first  to  apply'PRATT's  theories  to  a  large  project  and  he 
chose  to  measure  the  depth  of  compensation  from  the  physical  surface  of  the  earth  which  in¬ 
troduced  the  complication  that  the  surface  thus  determined  was  not  in  hydrostatic  equilibrium. 
Since,  however,  many  scientists  considered  it  to  be  a  necessary  condition  that  the  surface  of 
compensation  should  be  in  hydrostatic  equilibrium,  they  chose  to  reckon  the  depth  of  com¬ 
pensation  from  the  geoid.  it  is  interesting,  however,  that  to  a  large  degree  the  schools  ob¬ 
tained  the  same  result:  namely  that  the  depth  of  compensation  is  approximately  120  km.  At 
the  same  time  as  PRATT  G.  B.  AIRY  presented  a  somewhat  more  detailed  hypothesis  concer¬ 
ning  the  interior  of  the  earth.  AIRY  also  accepts  in  principle  the  isostatic  concept,  but  makes 
it  more  precise  by  giving  the  outer  body  of  the  earth  (sial),  a  density  of  2. 67,  which  floats 
on  the  inner  layer  (sima)  which  has  a  density  of  3.30.  The  continents  should  then  more  or 
less  float  like  icebergs  on  the  interior  masses  of  the  earth.  If  these  theories  were  correct  one 
could  find  in  them  a  good  basis  for  the  ideas  which  were  later  presented  by  the  German,  A.L. 
WEGENER.  According  to  him  the  South  American  and  African  continents  were  at  one  time  a 
single  continent,  but  have  since  then  separated.  This  hypothesis  of  the  movements  of  conti¬ 
nents  has  been  much  discussed,  but  no  successful  geodetic  measurements  have  yet  been  carried 
out  which  could  prove  or  disprove  the  hypothesis.  We  do  not  know  the  distance  at  the  present 
time  than  to  approximately  several  hundred  meters  and  the  movements  which  are  supposed  to 
occur  probably  are  essentially  less.  Therefore  it  will  be  some  time  before  this  theory  can  be 
checked,  if  in  fact  this  will  ever  be  possible.  The  actual  significance  of  isostasy  is  that  it 
offers  the  possibility  of  reduction  of  the  disturbance  effects  of  masses  lying  outside  the  geoid. 

A  disadvantage  of  the  isostatic  reduction  method  is  that  it  is  very  tedious,  but  in  spite  of  this 
isotatic  methods  have  come  into  wide  use.  For  example,  in  Sweden,  Rikets  Allmlnna  Kart- 
verk  has  performed  a  computation  of  a  suitable  geoid  for  Sweden  with  the  aid  of  isostatic 
reduction  methods. 

Among  other  theories  for  the  elimination  of  masses  outside  the  geoid  one  put  forward 
by  the  Russian  scientist  M.P.  RUDZKY  in  1905  can  be  cited.  He  sought  to  eliminate  the 
disturbing  masses  by  an  imaginary  transportation  of  masses  of  all  the  material  outside  the 
geoid.  In  this  original  manner  RUDZKY  obtained  a  model  earth  which  he  considered  to  be 
free  of  disturbing  effects.  A  number  of  other  theories  followed  and  the  state  had  almost 
been  reached  in  1950  where  each  geodesist  had  his  own  theory  for  the  computation  of  the 
geoid. 


The  cause  for  this  dilemma  can  perhaps  be  sought  for  along  various  lines.  First  and 
foremost  the  situation  was  such  that  the  various  theories  were  often  so  copplicated  that  it 
could  be  difficult  for  any  critics  to  find  the  correct  points  of  attack.  Furthermore  all  the 
theories  contained  fundamental  hypotheses  concerning  the  structure  of  the  interior  of  the 
earth  which  were  impossible  to  check  completely. 

It  was  therefore  a  scientific  sensation  of  the  first  order  when  at  the  General  Assembly 
for  Geodesy  and  Geophysics  in  Toronto  in  1957,  a  paper  written  *  by  the  Russian  scientist 
M.  MOLODENSKY  was  presented  in  which  it  was  shown  that  the  classical  geoid  is  not  re¬ 
quired  and  that  a  mathematically  correct  method,  free  of  hypotheses,  developed  in  the 
Russian  language,  had  already  been  available  for  several  years.  Furthermore,  this  was  gene¬ 
rally  applied  in  the  Soviet  Union.  The  situation  is  such  that  for  a  flat  lowland  the  previous 
theories  can  be  applied  when  the  effects  there  are  negligible.  MOLODENSKY  had  computed 
the  effect  of  disturbance  masses  on  the  plumb  line  for  a  mountain  4000  metres  in  height. 

A  correct  computation  showed  a  plumb  line  deflection  of  50",  while  a  computation  according 
to  classical  method  gives  a  result  of  15,4".  It  is  seen  that  the  older  method  does  not  even 
give  the  correct  order  of  magnitude. 

The  greatest  significance  of  MOLODENSKY' s  work  certainly  does  not  lie  in  the  purely 
practical  field.  Above  all,  it  is  through  the  mathematical  concepts  themselves  developed 
correctly  and  elegantly  that  assumptions  are  afforded  for  reasonable  analyses  of  all  the  prob¬ 
lems  which  are  connected  with  the  gravitational  field  of  our  earth. 
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The  mass  M  according  to  C i  also  defines  a  so-called  gravitational  potential 
W  according  to  the  relation 

w  =  f  w  J dv  (4) 

r  =  distance  of  the  volume  element 
from  the  actual  point 
Ouantities:  gal.  cm.  and  kilogal- meters. 


The  gravitational  potential  is  considered,  unless  otherwise  stated,  as  a  purely 
mathematical  definition.  According  to  (4)  the  gradient  to  the  gravitational  potential  is 


grad  W  * 
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-  -'Mi 

cos  (r,  x) 
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cos  (r,  y) 
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cos  (r,  z) 
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. 

(5) 


As  is  seen  from  equations  (2)  -  (5).  grad  W  constitutes  die  acceleration  which  is 
caused  by  the  actual  mass.  From  this  it  also  follows  that  the  potential  difference  between 
two  points  represents  the  work  which  is  required  to  move  a  unit  mass  from  one  point  to 
the  other.  Generally,  one  can  say  that  the  potential  according  to  equation  (4)  represents 
the  work  required  to  move  a  unit  mass  from  the  actual  point  to  infinity.  This  "absolute 
potential"  is  often  comparatively  difficult  to  obtain  and  in  general  one  must  accept  a  cer¬ 
tain  constant  error  in  the  determination  of  W  . 


GRAVITY 


A  body  at  rest  is  not  influenced  by  any  gravitation  other  than  that  defined  by  NEWTON. 
A  moving  body  is  also  affected  by  the  so-called  motional  gravitation.  For  the  earth 
this  additional  gravitation  is  generated  on  the  first  hand  by  the  earth's  rotation  about  its  own 
axis.  This  centrifugal  force  (g)  is  computed  as  follows 


where 


2 

a;  p 


w  =  angular  velocity  in  radians 

p  =  distance  between  the  actual  point 
to  the  axis  of  rotation. 


axis  thus  is 


P  2  2  2 

Q  -  to  x  dx  =  0.5  up 

0 


If  the  rotating  body  has  a  mass,  the  total  potential  is 

9  ...  .  .  2  2 


W  =  f  dV  +  0.5  oTp 


This  composed  potential  is  generally  called  the  gravity  potential.  The  gradient 
to  this  potential  evidently  constitutes  the  acceleration  caused  by  gravity. 
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where 

cos  (p,  z)  =  0  for  the  co-ordinate  system  in  which  the  z-axis  coincides  with  the 
axis  of  rotation. 


At  each  point  the  gradient  is  normal  to  the  equipotential  surface  which  includes  the 
point.  Furthermore  the  gradient  is  tangent  to  the  "plumb  line".  It  is  also  evident  that 
horn  a  rigorously  mathematical  point  of  view  the  "plumb  line"  cannot  be  considered  as  a 
straight  line.  However  in  most  practical  applications  this  should  be  justified. 


THE  CONCEPT  OF  THE  GRAVIMETRIC  POTENTIAL 


The  concept  of  potential  has  been  extensively  applied  in  electrical  theory  for  a  long 
time.  In  modern  geodesy  the  concept  of  potential  also  plays  a  significant  role.  To  explain 
the  geodetic  concept  of  potential  in  greater  detail  it  may  be  useful  to  begin  with  the  better 
known  definition  of  work.  If  two  points  A  and  B  are  given,  it  is  possible  to  define 
uniquely  the  work  required  to  transport  the  mass  m  from  the  lower  point  A  to  the 
higher  point  B  . 

B 

Work  A  -  B  *  ^  mgdh 
A 


where 

g  *  the  acceleration  of  gravity  along  the  chosen  path 
db  »  height  differences 


Here  the  potential  difference  is  identical  with  the  work  required  to  transport 
a  unit  mass  along  the  same  path.  Consequently  we  get  the  potential  difference 


B 

A  -  B  «  $  gdh  ■  Z  gdh  ■  g1dh1  ♦  gjdhj  ♦  ggdhg  ♦  g^  ♦  gjdhj  ♦  ggdhg  ♦  g^ 


The  potential  difference  is  the  only  mathematical  expression  free  of  objections  for 
"height  differences"  between  two  points.  For  practical  purposes  one  of  the  following  height 
concepts  is  often  used. 


Unreduced  height 


In  levelling  a  direct  measurement  of  the  quantities  dh  is  made.  For  simpler 
measurements  the  variations  In  gravity  can  be  disregarded  and  the  height  differences  be 
between  two  points  be  given  as 


It  is  to  be  noted  however,  that  such  a  computation  does  not  give  a  unique  result 
since  the  result  is  directly  dependent  on  the  path  chosen  between  the  two  given  points. 


Orthometric  height 

The  international  gravity  formula  reads 

y  =  978.049  (1  +  0.0052884  sin2  <p  -  0.0000059  sin2  2 <p)  cmsek'2 

From  a  comparison  of  the  gravity  values  at  die  equator  and  pole  we  get 

rp  *  rE  <’*005288> 

Consequently',  if  two  potential  differences  of  equal  magnitude  are  compared  we  obtain 

y  •  h  *  y  •  h„ 

'p  P  EE 


b*  =  -*h 

6  rE  p 


1.0053  hp 


In  a  corresponding  manner  a  theoretical  correction  can  be  applied  to  points  situated 
at  arbitrary  latitudes.  The  cotrection  is  called  the  orthometric  correction.  The 
correction  assumes  that  the  actual  gravity  agrees  with  the  theoretical  gravity. 

Dynamic  height 

The  potential  difference  between  two  points  A  and  B  is  defined  as 

A  %  ^ 

To  decide  "in  which  way  the  water  will  run"  orthometric  heights  cannot  be  used. 
After  division  of  the  potential  difference  by  a  suitable  mean  value  of  gravity  we  get  a 
height  difference 


This  height  difference  is  called  "the  dynamic  height  difference”.  Only  the  dynamic 
height  difference  give  the  so-called  "work  heights". 


THE  THEORETICAL  MODEL  EARTH 


For  the  geometric  definition  of  the  theoretical  model  earth  an  ellipsoid  (1924)  with 
die  following  data  is  used. 

Semi-major  axis:  a*  6  378  388.000  m. 

Semi-minor  axis:  b  =  6  356  911.946  m. 

a  -  b  1 

Flattening:  *  — —  =  a 

6  a  297 

-4 

Angular  velocity:  u  =  0.729  211  513  *  10  /sec.  (radians) 

For  the  geophysical  definition  of  die  model  earth  the  following  data  is  used. 

Gravity  on  the. 

ellipsoid's  surface:  y  =  978049  (1  +  0.0052884  sin  <p 

0  2 

-  0.0000059  sin  2  <p)  (milligal) 

where 

tp  •  latitude. 

\ 

Gravity  at  the  elevation  z 

dy  1  dy  2 

outside  the  ellipsoid:  y  ■  y  ♦  *r*  z  +  —  — *  z  ♦  ...  or 

r  z  o  8z  2  .  2 

dz 

y  ■  yQ  -  (308.78  -  0.45  sin2  <p)  z  +  0.0727  z 

2  .  z  (milligal  and  km) 

yz  ■  X0 - 1 —  ♦  ...  (f«  *he  sphere) 

Theoretical  height  z 

from  the  potential 

difference  (W  -  W) 
o 

WWW-W2  W-W3 

— -  ♦  — -  (0.000  157  854  -  0.000  001  034  sin2  <t)  ♦  — - 

y  V  y 

o  fo  fo 


z  * 


0.000000025 
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THE  EARTH’S  FORM  FROM  POTENTIALS  AND  GRAVITY 


For  the  computation  of  the  shape  of  the  earth  by  means  of  gravity  measurements,  it  is 
assumed  that  the  acceleration  of  gravity  is  known  at  all  points  of  the  earth’s  surface.  It  is 
further  assumed  that  the  latitudes  and  longitudes  of  all  ppints  of  measurement  are  known,  while 
the  elevation  of  the  points  above  the  international  reference  ellipsoid  are  unknown  and  sought 
quantities.  The  solution  to  the  problem  offers  various  difficulties  and  therefore,  in  general,  a 
direct  determination  has  not  been  carried  out,  but  instead  the  choice  has  been  to  determine  the 
potential  of  a  theoretical  auxilliary  model  which  is  assumed  to  have  the  same  surface  as  the 
actual  earth.  Such  a  solution  also  assumes  a  knowledge  of  the  potential  differences  between  all 
the  measured  points  concerned.  In  this  case  one  can  obtain  from  GREEN’s  second  theorem  two 
integral  equations  which  define  the  relationship  between  the  potentials  and  gravity  of  the  two 
models. 

For  the  actual  earth  we  obtain 


where 


W_ 


±  ((  (W  jl  .  i  ds  (i) 

2*  1  dn  r  dn  ' 


Wp  *  potential  at  the  actual  point 

W  =  potential  at  the  running  point 

r  =  distance  between  the  actual  point 

running  point 

n  ■  normal  to  the  surface  (outer) 

S  =  surface. 


and  the 


For  die  theoretical  model  earth  we  obtain 
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a- 


2ii  11  3n  r  dn 


where 

11^  *  theoretical  potential  at  the  actual  point 
U  =  theoretical  potential  at  the  running  point. 


(2) 


Neither  of  the  two  integral  equations  are  correct  for  the  case  represented  here  since 
GREEN's  theorem  requires  that  the  potentials  be  harmonic  functions,  e.g.  that  the  following 
LAPLACE  equations  should  be  satisfied 
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(8) 
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dx  dy  dz 


These  LAPLACE  equations  are  only  satisfied  for  a  body  which  does  not  rotate. 
However,  in  the  further  computations  the  differences  between  the  two  integral  equations 
are  considered  and  die  potential  differences'  are  certainly  harmonic  functions. 
Consequently  we  obtain 

W  -  U  *  T  =  r-  \\  [~T  *  “  (T^  ‘  f^>l  ds 

p  p  p  2n  •>  L  dn  r  dn  dn  J 


•  disturbance  potential  at  the  actual  point 
T  =  disturbance  potential  at  the  running  point. 

We  now  want  to  solve  this  integral  equation  with  respect  to  T  .  We  obtain 
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-  g  cos  (g,  n) 


y  cos  (y,  n) 


where  n  «  unit  vector  of  the  normal. 


For  an  approximately  spherical  surface  we  consider  the  following  permissible  approxi¬ 


mations 


37  “  '  8 


Here  g  is  die  gravity  on  the  surface  of  the  earth  while  y  It  the  corresponding 
quantity  for  the  model  earth.  For  the  computation  of  the  theoretical  gravity  values  we 
require  first  a  knowledge  of  the  gravity  field  of  the  model  earth.  From  the  geometric  point 
of  view  the  model  earth  is  assumed  to  be  a  rotational  ellipsoid  with  the  semi-major  axis  of 
6  378  388  m  and  a  flattening  of  1:297.  Gravity  on  the  surface  of  the  ellipsoid  is  assumed 
to  be  defined  according  to  the  International  Gravity  Formula  at  follows 

y  ■  978049  (1  ♦  0.0082884  sin2  ?  -  0.0000069  sin2  2?)  mgal  (10] 


^  ■  latitude. 


Gravity  at  the  elevation  z  above  th>:  ellipsoid  can  be  expressed  by  means  of  a  series 
development  of  the  following  type 


dy  1  d2)  2  1  d3y 
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dz  dz 
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For  the  theoretical  model  earth  we  get 

r 


=  y  -  (308.18  -  0.45  sin2  <p)  z  ♦  0.0121  z2  +  .. 
z  0 


or  for  the  sphere 

y , 


2y  z 
'0 


+  ... 


(Units:  milligal  and  km.) 


(R  *  radius  of  sphere) 


(11) 


(12) 


(12  a) 


Since  the  height  of  the  physical  surface  of  the  earth  above  the  ellipsoid  is  an  unknown 
quantity  it  is  not  possible  to  perform  a  direct  computation  of  the  theoretical  gravity  according 
to  (10).  However,  it  is  assumed  here  that  the  actual  potentials  W  are  known  for  every  point 
on  the  surface  of  the  earth.  Furthermore  it  is  assumed  that  the  theoretical  model  earth  is  an 


equipotentlal  surface  with  a  known  potential  U 


The  potential  differences  W  -  define 


a  height  difference  in  the  potential  field  of  the  model  earth  as  follows 
W  «  U  + 


‘Jk  •  3 

dz  *  2  2  z  *9  3  z  + 

dz  dz 


(13) 


A  known  potential  difference  W  -  Uq  can  be  evaluated  as  a  "theoretical  hei^tt" 


according  to  the  formula 


U  -  W  /u  -  w\8  /u  -  w\3 
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(14) 


For  the  model  earth,  we  obtain  the  corresponding" theoretical  height" 


U  -  W  /u  -  wV 


000  151  854  -  0.000  001  034  sin 


/u  -  w\8 

"'far) 


0.000000025 ♦ 
(15) 


Units:  U^  and  W  in  kilogal  and  meters;  y^  in  gai;  z  in  km. 


Here  it  is  assumed  that  the  z-axls  at  each  point  is  orthogonal  to  the  ellipsoid. 

Computation  of  the  corresponding  gravity  (y  )  can  now  be  performed  simply  according 

z 

to  equation  (12) 


dU 


Consequently  the  normal  derivative  *7—  can  be  denoted  by 

da 
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(16) 


For  the  sphere  the  corresponding  expression  can  be  denoted  as 
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(16  a) 


R  =  radius  of  the  earth. 

From  equations  (S)  and  (16)  we  obtain 
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Here  we  have  an  Integral  equation  with  T  as  the  unknown  quantity.  After  T  is 
computed  the  corresponding  height  difference  (Az)  can  be  determined  according  to  equation 
(IS).  The  final  height  above  the  ellipsoid  is  (z  +  Az). 

For  a  spherical  surface 
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Consequently  die  integral  equation  for  an  approximately  spherical  surface  is 
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(19) 


(19  a) 


This  integral  equations  can  be  solved  by  means  of  spherical  functions.  The  following 
parameter  is  introduced 


co 


T  *  E  E  a  Y  (tp,  X)  =  EE  a  Y 
n-0  |m|«1  nm 


(20) 


Thus  we  get 

EE  * Y  -  * Y  EE -srr, Y  Y) dS  =  i EE Y  Y <8-  >z> ds  <21> 

A  further  parameter  is  Introduced. 

b "  £  ss  £rr  -  v ds 

and  from  Integral  equation  (21)  we  obtain 
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The  relationship  between  a  and  b  can  now  be  determined 

l2n  a  =  b 
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or 
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Finally  we  obtain 
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2n  +  1 
2(n-  1) 


b  Y 


(2n  +  l)  Y 
2  (n  -  1 


<A 


2  Y 
2n  +  1 


(g  -  rz)  dS 


(24) 


T  =  $$  EE  — ~-r  Y  Y  (g  -  y  )  dS 
' 7  n  •  i  z 


(25) 


In  terms  of  LEGENDRE  polynomials  we  obtain 

T  =  iiSS  S(7fr,  Pn(C0SW)(g-  y2)dS  (26) 

Consequently  the  disturbance  potential  can  now  be  computed  horn 

T  -  ft  k  (g  -  Xz)  dS  (27) 

where 

k  »  jj-  (cotec  ^  -  6  tin  ^  ♦  l  -  cot  w  Ts  +  3  In  (tin2  ^  ♦  sin  ■*)”)  )  (28) 


This  it  the  well  known  STOKES’  formula. 


equipotential  surfaces 


A  surface  on  which  the  potential  is  everywhere  die  same  is  called  an  equipotential 
surface.  According  to  GREEN's  formula  the  following  is  valid  on  the  surface  of  a  body 


or 


a! 


W  ■  (W  -  -  •—)  dS 

p  2ir  >'  da  r  dn 

Thus  for  an  quipotential  surface  we  obtain 

WP  [’  Wt  cos  (r**  “> ds]  =  •  £  H  7  a?  dS 
WP  [’  +5T  T  cos  <*’  'n)  ds]  =  5T  Wf  005  <*•  *n) 


dS 


For  the  sphere 


cos  (r ,  n) 


—  and  we  get 


W 

P 


Finally,  if  g  is  also  constant  we  obtain 


From  the  geodetic  point  of  view  the  most  interesting  equipotential  surface  is  that 
equilibrium  surface  which  coincides  with  the  oceans  of  the  world.  In  geodesy  this  equi¬ 
potential  surface  is  called  the  geoid. 


MATHEMATICAL  PROPERTIES  OF  AN  EQ  UIPO  TE  N  T IA  L  SURFACE 


A  level  surface  is  characterized  by  the  fact  that  at  each  point  on  die  surface  the  potential 
is  of  the  same  magnitude. 

if  we  denote  W  as  the  scalar  expression  for  the  potential  we  have 

W  =  constant. 

in  an  arbitrary  point  Pq  the  potential  is  VI  ^  . 

If  Pq  is  made  the  origin  in  a  coordinate  system  we  get  the  following  expression  for 
the  potential  in  an  arbitrary  point  after  a  TAYLOR  series  development 
w  (x,  y.  z)  -  Wo  +  xwx  ♦  ywy  ♦  z ♦ 

♦  I  J  x2W  ♦  y2W  ♦  Z2W  +  2xy  W  +  2xz  W  ♦  tfi  W 

2  1  tx  '  yy  zz  xy  xz  yz 


♦  terms  of  higher  order 


where 


W 

x 

W 

y 

etc. 


partial  derivative  with  respect  to  x 

-  "  -  y 


VsTth  the  aid  of  matrices  this  expression  becomes 
W  (x.  y,  z)  -  Wo+  [x  y  zj 
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+  terms  of  higher  order. 

Here  the  potential  is  expressed  by  mean*  of  the  following  parameters 


1 :0 

Potential  scalar  - 

w 

0 

2:o 

Potential  vector  = 

W 

X 

W 

y 

*.] 

3:0 

Potential  matrix  = 

W 
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xy 
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W 

ZZ 

POTENTIAL  VECTORS  AND  PLUMB  LINES 


The  potential  vector  la  alio  denoted  by  the  gradient  to  the  level  turface.  For  this  the 
following  lymbols  are  used 


VW  ■  grad  W  « 


’aw" 
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X 

ax 

aw 
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y 

dy 

aw 

dz 

w 

z 

g 


From  this  we  obtain  the  relation 

g  •  g  =  W2  ♦  W*  ♦  W2  *  g2-  l  *=  g2 

It  ii  evident  that  the  gradient  it  tangent  to  the  plumb  line  at  the  actual  point.  We  orient 
the  coordinate  axes  10  that  die  z-axis  coincides  with  the  normal  to  the  level  surface  and  obtain 

W  ■  0,  W  *  0  and  W  =  a 
*  y  z  6 

The  radius  of  curvature  of  die  plumb  line  is  computed  from  die  principal  normals  as 
Sallows 


In  die  actual  case  we  get 
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From  this  we  obtain  for  the  normal  to  the  plumb  line 


and 


xx 


y* 


1 


where  c  =  the  derivative  of  the  acceleration  for  the  direction  <p  where  tg  ip  * 
°n 


POTENTIAL  MATRICES  AND  THE  LEVEL  SURFACE 


We  have  already  oriented  the  coordinate  system  so  that  the  z-axis  coincides  with  the 
normal  to  the  level  surface.  We  complement  this  with  a  reorientation  of  the  xy*plane  so 
that  die  rectangular  xy-terms  are  eliminated. 

For  this  we  determine  the  latent  roots  to  the  system 


which  gives 
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The  potential  equation  can  now  be  written 

W  (x  y  z)  -  W#tZW2  ♦  Y8  X8  ♦  Z2Wa)  ♦  ...  ♦ 


♦  terms  of  hitter  order 


The  quadratic  potential  term  can  now  be  easily  ev  aimed  in  terms  of  height  differences 
after  division  by  g  . 


dZ., 


X% 


dZ. 


Y\ 

a 

2g 


where 


B  -8- 

R, 


r  A 
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R  and  R  are  here  the  two  principal  radii  of  curvature  to  the  level  surface. 

'  & 

(The  X-  and  Y-axes  are  neceuarily  orthogonal  in  this  system.) 

The  trace  of  the  complete  potential  matrix  is  known  as  LAPLACE's  equation 


SpW  =  W  ♦  W  +  W  «  AW 
xx  yy  zz 

where  Sp  W  ■  trace  of  W  . 

For  NEWTON's  gravitation  in  empty  space 

Sp  W  -  0  is  valid. 

A  potential  which  satisfies  this  condition  is  said  to  be  harmonic. 
When  the  actual  body  rotates  we  get 

2 

Sp  W  =  2w  w  «  angular  velocity 

and  where  the  density  of  the  mass  *  9 

2 

Sp  W  *  2w  -  4»  f  9 


f  ■  NEWTON’s  gravitational  constant. 


SIMPLE  LAYER  POTENTIAL 


The  potential  is  a  unique  and  continuous  scalar  quantity.  Ignoring  the  choice  of  units, 
the  potential  according  to  NEWTON  can  be  denoted  as  follows 


9 

r 


W  =  % 
V 


3 

density  of  mass  (kg/m  );  V  =  volume 
distance  between  volume  element  and  the  actual  point 


In  the  interior  of  the  mass  with  continuous  mass  distribution 

o 

7  W  =  -  4a  9  (POISSON's  equation) 
is  valid  and  in  empty  space 

V2  W  =  0  (LAPLACE's  equation) 


In  die  latter  case  the  function  is  said  to  be  harmonic. 

While  the  above  "volume  potential"  is  connected  directly  to  known  physical  rela¬ 
tions,  parallel  to  this  we  use,  for  example,  the  potential  from  a  hypothetical  layer  of  mass  on 
die  surface. 


W  *  M  ~  dS  V2  W  =  -  4 a  o 

2 

o  =  layer  density  (kg/m  );  S  *  surface 

r  -  distance  between  the  surface  element  and  the  actual  point. 

In  this  case  we  employ  a  so  called  single  layer  potential.  The  procedure  lacks 
physical  background  in  geodesy  but  nevertheless  constitutes  a  valuable  mathematical  auxilliary 
tool. 

While  the  gradients  to  the  volume  potential  are  both  unique  and  continuous,  the  single 
layer  potential  derivatives  exhibit  a  "pimp"  (discontinuity)  at  the  surface.  This  can  be  shown 
in  the  following  manner. 


Consider  an  infinitesimal  cylinder  which  cuts  through  the  surface  S  and  is  parallel  to 
the  surface  normal 


If  we  start  from  POISSON's  equation  we  obtain,  after  integration 

^  V^W  dV  =  *  ^  4*  o  dV 
V  V 


From  GAUSS'  theorem  and  GAUSS’  law  we  further  obtain 


\  V  W  dV  ■  $  W  dS  =  -  4*  ^  o  dS 

VS  S 

If  the  gradients  are  computed  at  right  angles  to  the  surface,  they  are  parallel  to  the 
vector  dS  and  we  get 

aw  aw. 


From  this  we  get 


v6n  3n. 

y  1 


aw  aw 


)  *  -  4t  o 


dn 


an. 


4*  o 


To  deduce  an  expression  for  the  continuous  part  of  the  derivatives  we  differentiate  in 
the  direction  of  the  tangent 


aw 

dt 


8r 

at 


dS 


When  the  continuous  part  is  known  we  obtain  the  derivatives  in  the  direction  of  the 
normal  after  addition  of  the  discontinuities  above 


For  the  computation  of  the  derivative  in  an  arbitrary  direction  we  perform  a  vectorial 
addition  according  to  the  figure  below 


aw  aw  aw  ,  % 

—  «  —  coife  «n)  ♦—<»(«».  t) 


aw  /  a  rdr  dr  ~t 

7^  *  *  >  7  LTn  601  (n*  m)  *  Tt  001  (m*  ‘>J  «* '  ^o®0* m) 

S  r 


where 


—  =  cos  (r,  n) 

on 


dr 

3t 


cos  (rp  t) 


From  these  we  finally  obtain 


where 


(r,  m) 

(n,  m) 
a 

o 


angle  between  the  vector  m  and  the  vector  from  P^  to  the 
running  point 

angle  between  the  surface  normal  and  the  vector  m 
layer  density  at  the  point  . 
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